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ABSTRACT 

We report the discovery of the closest known black hole candidate as a binary companion to V723 Mon. V723 Mon is a 
nearby (d ~ 460pc), bright (V = 8.3 mag), evolved (Tem,sia ~ 4440 K, and Leian =~ 173 Lo) red giant in a high mass 
function, f (M) = 1.72 + 0.01 Mo, nearly circular binary (P = 59.9 d, e =~ 0). V723 Mon is a known variable star, previously 
classified as an eclipsing binary, but its All-Sky Automated Survey (ASAS), Kilodegree Extremely Little Telescope (KELT), 
and Transiting Exoplanet Survey Satellite (TESS) light curves are those of a nearly edge-on ellipsoidal variable. Detailed models 
of the light curves constrained by the period, radial velocities and stellar temperature give an inclination of SEE n a mass 
ratio of q = 0.33 + 0.02, a companion mass of Mcomp = 3.04 + 0.06 Mo, a stellar radius of Reian = 24.9 € 0.7 Ro, and a 
giant mass of Mgiant = 1.00 + 0.07 Mo. We identify a likely non-stellar, diffuse veiling component with contributions in the 
B and V-band of ~63% and ~24%, respectively. The SED and the absence of continuum eclipses imply that the companion 
mass must be dominated by a compact object. We do observe eclipses of the Balmer lines when the dark companion passes 
behind the giant, but their velocity spreads are low compared to observed accretion disks. The X-ray luminosity of the system 
is Lx = 7.6 x 10? ergs s^!, corresponding to L/Leaa~10~°. The simplest explanation for the massive companion is a single 
compact object, most likely a black hole in the *mass gap". 


Key words: stars: black holes — (stars:) binaries: spectroscopic — stars: individual: V723 Mon 


* E-mail: jayasinghearachchilage. 1 @osu.edu 


© 2021 The Authors 


2  T.Jayasinghe et al. 


1 INTRODUCTION 


The discovery and characterization of neutron stars and black holes in 
the Milky Way is crucial for understanding core-collapse supernovae 
and massive stars. This is inherently challenging, partly because iso- 
lated black holes are electromagnetically dark and partly because 
compact object progenitors (OB stars) are rare. To date, most mass 
measurements for neutron stars and black holes come from pulsar and 
accreting binary systems selected from radio, X-ray, and gamma-ray 
surveys (see, for e.g., Champion et al. 2008; Liu et al. 2006; Ozel 
et al. 2010; Farr et al. 2011), and from the LIGO/Virgo detections of 
merging systems (see, for e.g., Abbott et al. 2016, 2017). Interacting 
and merging systems are however a biased sampling of compact ob- 
jects. A more complete census is needed to constrain their formation 
pathways. 

One important component of such a census is to identify non- 
interacting compact objects in binaries around luminous compan- 
ions. By their very nature, interacting black holes only sample a 
narrow range of binary configurations, and almost the entire pa- 
rameter space of binaries with black holes that are non-interacting 
remains unexplored. Interacting compact object binaries are only 
active for relatively short periods of time, so most systems are quies- 
cent or non-interacting. The discovery and characterization of these 
non-interacting black holes are important for understanding the birth 
mass distribution of black holes and their formation mechanisms. 

With advances in time-domain astronomy and precision Gaia as- 
trometry (Gaia Collaboration et al. 2018), a significant number of 
these systems should be discoverable. For example, Breivik et al. 
(2017) estimated that ^10? — 104 non-interacting black holes are 
detectable using astrometry from Gaia. Similarly, Shao & Li (2019) 
used binary population synthesis models to estimate that there are 
—10? detached non-interacting black holes in the Milky Way, with 
10? of these systems having luminous companions that are brighter 
than G~20 mag. 

Thompson et al. (2019) recently discovered the first low-mass 
(MBH = Leu Mo) non-interacting black hole candidate in the 
field. It is in a circular orbit with Po, ~ 83d around a spotted 
giant star. Other non-interacting BH candidates have been discov- 
ered in globular clusters: one by Giesers et al. (2018) in NGC 3201 
(minimum black hole mass Mg = 4.36 + 0.41 Mo), and two by 
Giesers et al. (2019) in NGC 3201 (Mgy sin(i) = 7.68 + 0.50Mo 
and Mpysin(i) = 4.4 + 2.8 Mo). While obviously interesting in 
their own right, these globular cluster systems likely have formation 
mechanisms that are very different from those of field black hole 
binaries. 

Other claims for non-interacting BH systems have been ruled out. 
For example, LB-1, which was initially thought to host an extremely 
massive stellar black hole (Mgg = 68*1! Mo, Liu et al. 2019), 
was later found to have a much less massive companion that was 
not necessarily a compact object (see, for e.g., Shenar et al. 2020; 
Irrgang et al. 2020; Abdul-Masih et al. 2020; El-Badry & Quataert 
2020b). Similarly, the naked-eye star HR 6819 was claimed to be a 
triple system with a detached black hole with Mgy = 6.3 + 0.7 Mo 
(Rivinius et al. 2020), but was later argued to contain a binary system 
with a rapidly rotating Be star and a slowly rotating B star (El-Badry 
& Quataert 2020a; Bodensteiner et al. 2020). 

Here we discuss our discovery that the bright red giant V723 Mon 
has a dark, massive companion that is a good candidate for the closest 
known black hole. We discuss the current classification of this system 
in Section 2, and describe the archival data and new observations used 
in our analysis in Section 3. In Section 4, we analyze photometric 
and spectroscopic observations to derive the parameters of the binary 
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system and the red giant secondary. In Section 5, we discuss the 
nature of the dark companion. We present a summary of our results 
in Section 6. 


2 V723 MON 


V723 Mon (HD 45762, SAO 133321, TIC 43077836) is a luminous 
(my = 8.3) red-giant! in the Monoceros constellation with J2000 
coordinates (a, 6) = (97.269410?, —5.572286?). It was classified as 
a likely long period variable in the General Catalogue of Variable 
Stars (GCVS; Kazarovets et al. 1999) after it was identified as a vari- 
able source in the Hipparcos catalogue with period P = 29.97 d (ESA 
1997). Subsequently, the All-Sky Automated Survey (ASAS) (Poj- 
manski 1997, 2002) classified it as a contact/semi-detached binary 
with P = 59.87 d. The Variable Star Index (VSX; Watson et al. 2006) 
presently lists it as an eclipsing binary of the G-Lyrae type (EB) with 
P = 59.93d. 

V723 Mon has a well determined spectroscopic orbit and is in- 
cluded in the the Ninth Catalogue of Spectroscopic Binary Orbits 
(Pourbaix et al. 2004). In particular, Griffin (2010) identified V723 
Mon as a single-lined spectroscopic binary (SB1) with a nearly circu- 
lar P~60 d orbit. Strassmeier et al. (2012) (hereafter S12) refined the 
orbit to Porb = 59.9363 + 0.0016 d and e9,52:0.0150 + 0.0009 using 
a large number of high precision RV measurements from STELLA 
(Weber et al. 2008). S12 also argue that this is the outer orbit of 
a triple system, and that the more massive (inner) component is an 
SBI consisting of another giant star and an unseen companion with a 
period Pinner ~ Pouter/3. Griffin (2014) (hereafter G14) was unable 
to find a spectral feature indicative of a second companion in their 
cross-correlation functions. G14 discusses several peculiarities in the 
S12 RV solution. In particular, the radial velocity curve associated 
with the second component is unusual in structure compared to any 
other system characterized by $12 and a triple system with this period 
ratio would almost certainly be dynamically unstable. 

The most striking feature of the well-measured 60 day RV curve 
is its large mass function of 


Pob K? (L - e?)3/2 uH ME nis sin? i 
(27G) (Mgiant + Mcomp)? 


given Porb = 59.9363 d, e = 0.015 and K = 65.45 km s^! from S12. 
If the observed giant has a mass of Mgiant ~ 1 Mo, the mass function 
implies a massive companion with a minimum mass of Mcomp ^ 
3 Mo. Since the observed light is clearly dominated by the giant and 
the companion has to be both much less luminous and significantly 
more massive than the giant, V723 Mon is a prime candidate for 
a non-interacting, compact object binary. This realization led us to 
investigate V723 Mon in detail as part of a larger project to identify 
non-interacting, compact object binaries. 


f(M)- = 1.73 Mo, (1) 


1 V723 Mon has previously been assigned a spectral type of GO II (Houk & 
Swift 2000), however, in this work we find that it is more consistent with a 
KO0/K1 III spectral type. 


3 OBSERVATIONS 


Here we present observations, both archival and newly acquired, that 
are used in our analysis of V723 Mon. 


3.1 Distance, Kinematics and Extinction 


In Gaia EDR3 (Gaia Collaboration et al. 2020), V723 Mon is 
source. id-3104145904761393408. Its EDR3 parallax of WEDR3 = 
2.175 + 0.033 mas implies a distance of d = 1/w = 460 +7 pc. 
This is little changed from its Gaia DR2 parallax of «pg» = 
1.9130 + 0.0563 mas. At these distances, there is little difference 
between d = 1/« = 519 + 15 pc (for DR2) and the more careful 
estimate of d = 515.64136 pc by Bailer-Jones et al. (2018). The 
astrometric solution has significant excess noise of A — 0.22 mas, 
which is not surprising given that the motion of the giant should be 
— 0.8 mas. However, its renormalized unit weight error (RUWE) of 
1.39, while larger than unity, is not indicative of problems in the 
parallax. We adopt a distance of d = 460 pc for the remainder of the 
paper. The distance uncertainties are unimportant for our analysis. 

V723 Mon has Galactic coordinates (1, b) = (215.372°, —7.502?), 
close to the Galactic disk, but away from the Galactic center. At 
the EDR3 distance, V723 Mon is ~32 pc below the midplane. Its 
proper motion in EDR3 is us = —1.347 x 0.032 mas yr !, and 
us = 16.140 + 0.031 mas yr~!. Combining this with the systemic 
radial velocity from $4.2 and the definition of the local standard 
of rest (LSR) from Schónrich et al. (2010), the 3D space motion 
of V723 Mon relative to the LSR as (U, V,W)Lsr = (-10.8 + 
4.3, 36.9 + 5.3, 20.1 € 3.4) km s using BANYAN (Gagné et al. 2018) 
for the calculations. We calculated the thin disk, thick disk and halo 
membership probabilities based on the UVW velocities following 
Ramirez et al. (2013) to obtain P(thin) ~ 97%, P(thick) ~ 3% 
and P(halo) ~ 0%, respectively. This suggests that this system is a 
kinematically normal object in the thin disk. 

Gaia DR2 (Gaia Collaboration et al. 2018) also reports a lumi- 
nosity Lgaia = 165.5 + 6.5 Lo, temperature Ter Gaia = 4690*20 K, 
and radius RGaia = 183 Ro, for the star that are consistent with 
an evolved red giant. While Gaia DR2 does not report a value for 
the reddening towards V723 Mon, Gontcharov & Mosenkov (2017) 
reports E(B — V) ~ 0.10. The three-dimensional dust maps of Green 
et al. (2019) give E(B — V) = 0.10 + 0.04 at the Gaia distance, 
consistent with this estimate. 


32 Light Curves 


We analyzed well-sampled light curves from the All-Sky Auto- 
mated Survey (ASAS) and the Kilodegree Extremely Little Tele- 
scope (KELT), a densely sampled but phase-incomplete light curve 
from the Transiting Exoplanet Survey Telescope (TESS), BVR.I. 
light curves from the Remote Observatory Atacama Desert (ROAD) 
and a sparse ultraviolet (UV) light curve from the Neil Gehrels Swift 
Observatory. 

ASAS (Pojmanski 1997, 2002) obtained a V-band light curve of 
V723 Mon spanning from November 2000 to December 2009 (~54 
complete orbits). We selected 580 epochs with GRADE-A or GRADE-B 
for our analysis. V723 Mon clearly varies in the ASAS light curve, 
with two equal maxima but two unequal minima, when phased with 
the orbital period from S12. We determined the photometric period 
using the Period84 software (Lenz & Breger 2005). The dominant 
ASAS period of Pasas~29.9674 + 0.0138 d corresponds to Porb/2. 
Once the ASAS light curve was whitened using a sinusoid, we find 


V723 Mon 3 


an orbital period of 
Porb,asas = 59.9863 + 0.0551 d, (2) 


which agrees well with the spectroscopic periods from S12 and G14. 
Unfortunately, V723 Mon is saturated in the Automated Survey for 
SuperNovae (ASAS-SN; Shappee et al. 2014; Kochanek et al. 2017) 
images, and we could not use it to extend the time span of the V-band 
data. 

The KELT (Pepper et al. 2007) light curve contains 1297 epochs 
which we retrieved from the Exoplanet Archive”. The KELT Rx filter 
can be considered as a very broad Johnson R-band filter (Siverd et al. 
2012). However, there can be significant color corrections compared 
to a standard Johnson R-band filter for very blue and very red stars 
(Pepper et al. 2007; Siverd et al. 2012). KELT observations were 
made between September 2010 and February 2015 (~26 complete 
orbits). The dominant period in the KELT data (Pkgj 7229.9682 + 
0.0279 d) again corresponds to Porb/2. We find an orbital period of 


Porb,KELT = 60.0428 + 0.1121 d. (3) 


The difference between the ASAS and KELT photometric period 
estimates is not statistically significant. 

V723 Mon (TIC 43077836) was observed by TESS (Ricker et al. 
2015) in Sector 6, and the 27 days of observations correspond to 
[0.46,0.82] in orbital phase where the phase of the RV maximum 
is 0.75. V723 Mon was also observed in Sector 33, with the ob- 
servations spanning [0.94,0.36] in orbital phase. We analyzed the 
TESS data using the adaptation of the ASAS-SN image subtraction 
pipeline for analyzing TESS full-frame images (FFIs) described in 
Vallely et al. (2020). While this pipeline produces precise differential 
flux light curves, the large pixel scale of TESS makes it difficult to 
obtain reliable measurements of the reference flux of a given source. 
We normalized the light curve to have the reference TESS-band 
magnitude of 7.26 in the TESS Input Catalog (Stassun et al. 2019). 
Conveniently, the mean of the Sector 6 observations is approximately 
the mean for a full orbital cycle (see Figure 5). We use a zero point 
of 20.44 electrons (TESS Instrument Handbook?). The light curve 
does not include epochs where the observations were compromised 
by the spacecraft's momentum dump maneuvers. 

We obtained BVR.I. light curves at the Remote Observatory 
Atacama Desert (ROAD; Hambsch 2012). All observations were ac- 
quired through Astrodon Photometric filters with an Orion Optics, 
UK Optimized Dall Kirkham 406/6.8 telescope and a FLI 16803 
CCD camera. Twilight sky-flat images were used for flatfield cor- 
rections. Reductions were performed with the MAXIM DL program* 
and the photometry was carried out using the LesvePhotometry 
program." 

We obtained Swift UVOT (Roming et al. 2005) images in the 
UVM2 (2246 A) band (Poole et al. 2008) through the Target of 
Opportunity (ToO) program (Target ID number 13777). We only 
obtained images in the UVM2 band because the Swift UVW1 and 
UVW2 filters have significant red leaks that make them unusable 
in the presence of the emission from the cool giant, and the star is 
too bright to obtain images in the optical UVOT filters. Each epoch 
of UVOT data includes multiple observations per filter, which we 
combined using the uvotimsum package. We then used uvotsource 


2 https: //exoplanetarchive.ipac.caltech.edu/ 

3 https://archive.stsci.edu/files/live/sites/mast/files/ 
home/missions-and-data/active-missions/tess/_documents/ 
TESS_Instrument_Handbook_v0.1.pdf 

^ http: //www.cyanogen. com 

5 http: //www.dppobservatory.net/ 


MNRAS 000, 1—27 (2021) 


4 T. Jayasinghe et al. 


to extract source counts using a 5"0 radius aperture centered on the 
star. We measured the background counts using a source-free region 
with radius of 40"0 and used the most recent calibrations (Poole 
et al. 2008; Breeveld et al. 2010) and taking into account the recently 
released update to the sensitivity correction for the Swift UV filters®. 
The Swift UV M2 observations are summarized in Table 6. We report 
the median and standard deviation of these UV M2 observations in 
Table 1. 


3.3 Radial Velocities 


We used two sets of radial velocity (RV) measurements. The first 
set, from Griffin (2014), was obtained between December 2008 and 
November 2013 as part of the Cambridge Observatory Radial Ve- 
locity Program and span 1805 days. The median RV error for this 
dataset is 0.80 km s^! . These 41 RV epochs were retrieved from the 
Ninth Catalogue of Spectroscopic Binary Orbits (SB9; Pourbaix et al. 
2004)," converting the reported epochs to Barycentric Julian Dates 
(BJD) on the TDB system (Eastman et al. 2010) using barycorrpy 
(Kanodia & Wright 2018). 

The second set of RV data consists of 100 epochs obtained by 
$12 with the high resolution (R ~ 55000) STELLA spectrograph. 
STELLA spectra have a wavelength coverage of 390-880 nm and 
a spectral resolution of 0.12 A at 650 nm. Spectra were obtained 
between November 2006 and April 2010, spanning a baseline of 1213 
days. The spectra were reduced following the standard procedures 
described in Strassmeier et al. (2012) and Weber et al. (2008). Of the 
100 spectra, 75 had S/N » 30 near 650 nm. There were 87 epochs 
with good RV measurements for the giant and the median RV error 
is ~0.19kms~!. The STELLA RV measurements are listed in Table 
T. 


3.4 Additional Spectra 


To better understand the V723 Mon system, and to test for possi- 
ble systematic errors, we obtained a number of additional high and 
medium resolution spectra These observations are summarized in 
Table 8. Using the HIRES instrument (Vogt et al. 1994) on Keck 
I, we obtained 7 spectra with R ~ 60000 between Oct 20 2020 
and Dec 26 2020 using the standard California Planet Search (CPS) 
setup (Howard et al. 2010). The exposure times ranged from 35 to 
60 seconds. We also obtained a very high resolution (R « 220000) 
spectrum on 29 Nov. 2020 using the Potsdam Echelle Polarimet- 
ric and Spectroscopic Instrument (PEPSI; Strassmeier et al. 2015) 
on the Large Binocular Telescope. We used the 100 um fiber and 
6 cross-dispersers (CD). The data were processed as described in 
Strassmeier et al. (2018). The total integration time was 90 minutes, 
and the combined spectrum covers the entire wavelength range acces- 
sible to PEPSI (3840 — 9070 À). The spectrum has a signal-to-noise 
ratio (S/N) of 70 in the wavelength range 4265 — 4800 A (CD2) and 
S/N= 260 in the range 7419 — 9070 A (CD6). 

Using the medium resolution (R ~ 2000) Multi-Object Double 
Spectrographs mounted on the twin 8.4m Large Binocular Telescope 
(MODS1 and MODS2; Pogge et al. 2010), we obtained a series of 
spectra from Nov 18 to Nov 22 2020 as the dark companion moved 
into eclipse behind the observed giant. Exposure times were typically 
15 seconds, but this was extended to 24 seconds for the Nov 22 obser- 
vation to compensate for clouds. We used a 2’’4 wide slit to ensure 


6 https://www.swift.ac.uk/analysis/uvot/index.php 
7 https://sb9.astro.ulb.ac.be/DisplayFull.cgi?3936+1 
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that all the light was captured. We reduced these observations using 
a standard combination of the modsccdred? pyrHon package, and 
the modsidl pipeline?. The blue and red channels of both spectro- 
graphs were reduced independently, and the final MODS spectrum 
for each night was obtained by averaging the MODS1 and MODS2 
spectra. The best weather conditions during this observing run oc- 
curred on Nov 20. The HIRES, PEPSI and MODS observations are 
summarized in Table 8. 


3.5 X-ray data 


We analyzed X-ray observations from the Swift X-Ray Telescope 
(XRT; Burrows et al. 2005) and XMM-Newton (Jansen et al. 2001). 
The XRT data were taken simultaneously with the UV M2 obser- 
vations, with individual exposure times of 160 to 1015 seconds. 
Additionally, two longer ~5 ks XRT exposures were taken on 2021- 
01-21 (9 = 0.5) and 2021-02-20 (9$ = 0). In total, V723 Mon was 
observed with Swift XRT for 19865 seconds. All XRT observations 
were reprocessed using the Swift xRTPIPELINE version 0.13.2 and 
standard filter and screening criteria! and the most up-to-date cali- 
bration files. To increase the signal to noise of our observations, we 
combined all cleaned individual XRT observations using XSELECT 
version 2.4g. To place constraints on the presence of X-ray emission, 
we used a source region with a radius of 30 arcsec centered on the 
position of V723 Mon and a source-free background region with a 
radius of 150 arcsec located at RA = 06:28:53.1, Dec =—05:38:49.6 
(J2000). 

We also retrieved archival XMM-Newton data obtained during 
a ~10 ks observation of the nearby ultraluminous infrared galaxy 
IRAS 06269-0543 (Observation ID 0153100601; PI: N. Anabuki). 
However, V723 Mon is ~12’ off-axis in these observations, resulting 
in a non-optimal PSF with a 90% enclosed energy radius of ~1’. 
We reduced the data using the XMM-Newton Science System (SAS) 
Version 15.0.0 (Gabriel et al. 2004). We removed time intervals 
with proton flares or high background after identifying them by 
producing count-rate histograms using events with an energy between 
10-12 keV. For the data reduction, we used the standard screening 
procedures and the FLAGS recommended in the current SAS analysis 
threads!! and XMM-Newton Users Handbook!”. 


4 RESULTS 


Here we present our analyses of the observations described in §3. 
In §4.1, we characterize the red giant using its spectral energy dis- 
tribution (SED) and spectra. In §4.2, we fit Keplerian models to the 
radial velocities and derive a spectroscopic orbit for V723 Mon. In 
§4.3, we model the ellipsoidal variations of the red giant using multi- 
band light curves and the binary modeling tool PHOEBE to derive the 
masses of the red giant and the dark companion. We also derive lim- 
its on companion eclipse depths. In §4.4, we characterize the veiling 
component in this system. In §4.5, we place constraints on luminous 


8 http://www.astronomy.ohio-state.edu/MODS/Software/ 
modsCCDRed/ 

9 http://www.astronomy.ohio-state.edu/MODS/Software/ 
modsIDL/ 

10 http://swift.gsfc.nasa.gov/analysis/xrt swguide vl 2. 
pdf 

11 https://www.cosmos.esa.int/web/xmm-newton/sas-threads 
12 https://xmm-tools.cosmos.esa.int/external/xmm user. 
support/documentation/uhb/ 


stellar companions. In $4.6, we characterize the Balmer emission 
lines and their variability. In $4.7, we discuss the X-ray observations. 


4.1 Properties of the Red Giant Secondary 


We characterized the red giant using both fits to its overall SED and 
analyses of the available spectra. For the SED we used photometry 
from APASS DR10 (Henden et al. 2018), SkyMapper DR2 (Onken 
et al. 2019), 2MASS (Skrutskie et al. 2006) and AIIWISE (Wright 
et al. 2010). We used the Swift UVM2 photometry only as an upper 
limit ($3.5). The compilation of the multi-band photometry used in 
these fits are given in Table 1. 

We fit the SED of V723 Mon using DUSTY (Ivezic & Elitzur 
1997; Elitzur & Ivezić 2001) inside a MCMC wrapper (Adams & 
Kochanek 2015). We assumed foreground extinction due to Ry = 
3.1 dust (Cardelli et al. 1989) and used the ATLASO Castelli & 
Kurucz (2003) model atmospheres for the star. We assume that the 
source is at the Gaia EDR3 distance and used minimum luminosity 
uncertainties of 19% for each band to obtain x Ndof = 1. The 
expanded uncertainties needed to reach y2/Nao f = 1 with respect 
to the reported photometric errors in each measurement are likely 
driven by using single epoch photometry for a variable source plus 
any systematic problems in the models and photometry. We used 
a weak Te giant = 4500 + 500 K prior on the temperature and a 
prior of E(B — V) = 0.10 + 0.04 on the extinction from Green 
et al. (2019). The SED fit yields Tog giant ~ 4440 + 90 K, Lgiant = 
173 + 8 Lo, Rgiant = 22.2 € 0.8 Ro and E(B — V) = 0.085 x 0.034 
mag. Figure 1 shows the SED and the best fitting model. The SED 
well-constrains the temperature and is consistent with the extinction 
estimates. If we force the model to have smaller stellar radii, the 
goodness of fit worsens rapidly, from y? = 13.6 for the best fit to 
20.3 for Rgiant = 20Ro and 40.9 for Rgiant = 18Ro. Decreasing 
the radius forces the star to become hotter (4800 K for 18Rọ) to 
fit the SED at long wavelengths and more extincted to fit it at short 
wavelengths. Constraints on a stellar companion to the giant from 
the SED are discussed in Section $4.5. 

The equivalent width (EW) of the Nar D doublet provides an 
independent estimate of the Galactic reddening towards V723 Mon 
(see for e.g., Poznanski et al. 2012). Using the very high resolution 
PEPSI spectrum, we measure EW(D4) = 125.1 + 12.0 mA and 
EW(D») = 85.7 + 7.5 mA. From the equivalent width—E(B — V) 
calibration in Poznanski et al. (2012), we find that E(B—V) = 0.025x 
0.005 mag. This is consistent with the low foreground extinction 
that was derived from the SED models and the Green et al. (2019) 
extinction maps. 

All high resolution spectra indicate that the giant is rapidly rotat- 
ing. Strassmeier et al. (2012) derived a projected rotational velocity 
of Vrot Sin? = 16 + 2kms ^! using the SES spectra. Griffin (2014) 
found a similar average vrot sini ~ 15 km s^! but noted that the values 
seemed to depend on the orbital phase and ranged from vrot sini ~ 
10 km s^! to Vrot sini ~ 20km s^. We also find that Vrot Sin į varies 
from ^17.6 km s^! to ~21.7 kms! in the SES spectra. The HIRES 
CPS pipeline (Petigura 2015) reports vrot sini = 21.6 + 1.0 km s7! 
and we found vrot sini = 17.9 + 0.4kms™! from the PEPSI spec- 
trum using iSpec (Blanco-Cuaresma et al. 2014; Blanco-Cuaresma 
2019). Assuming that the rotation of the giant is tidally synchronized 
with the orbit, the SES, PEPSI and HIRES measurements yield stel- 
lar radii of R sin i = 19.0 + 2.4 Rọ, Rsin i = 21.2 + 0.5 Ro and 
R sin i = 25.6 + 1.2 Ro, respectively, consistent both with estimates 
from the SED and the radius derived for the giant from the PHOEBE 
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Figure 1. The best-fitting, extinction-corrected SED model for V723 Mon 
without considering veiling (see §4.4). The Swift UVM2 detection (arrow; 
see Table 1) is treated as an upper limit and the error bars are expanded to 
give y7/ Naof =~ 1. The SEDs for a main sequence companion of mass 
1.29 Mo, an equal mass binary consisting of two main sequence stars each 
with 0.99 Mo, and a binary with companion masses of 1.25 Mo and 0.77 Mo 
are also shown (see $4.5). These SEDs closely overlap and are hard to distin- 
guish. As shown by the SED of a single 2.5 Mo main sequence star (dashed 
line), hotter stars are easily ruled out. 


models ($4.3). The match to the estimate of the radius from the SED 
indicates that sini ~ 1 independent of the PHOEBE models. 

To derive the surface temperature (Teg), surface gravity (log(g)), 
and metallicity ([Fe/H]) of the giant, we use the spectral synthesis 
codes FASMA (Tsantaki et al. 2020; Tsantaki et al. 2018) and iSpec. 
FASMA generates synthetic spectra based on the ATLAS-APOGEE 
atmospheres (Kurucz 1993; Mészáros et al. 2012) with MOOG (Sneden 
1973) and outputs the best-fit stellar parameters following a y? min- 
imization process. iSpec carries out a similar minimization process 
with synthetic spectra generated by SPECTRUM (Gray & Corbally 
1994) and MARCS model atmospheres (Gustafsson et al. 2008). 
The line lists used in this process span the wavelength range from 
480 nm to 680 nm. Since we are confident the companion should un- 
dergo eclipses (see $4.5), we fit the SES spectrum near conjunction 
(0.48 S ¢ < 0.52) when any companion would be eclipsed by the 
giant. For the detailed fits, we fix the rotational velocity to the value 
of vrot sin į = 19.1 + 0.4km s^! found by iSpec for this spectrum. 

For the FASMA fits we initially keep the macroturbulent broadening 
(Vmac) and the microturbulence (Vmicro) fixed at Vmac = 5kms7! 
and Vic = 2kms~!, but then allow them to be optimized once 
we have a reasonable fit. This fits yield Teff giant = 4480 + 50 K, 
log(g) = 1.7 + 0.2, and [Fe/H] = —1.1 + 0.1. In the iSpec fits, 
Vmic Was kept as a free parameter and we obtain similar results with 
Teff giant = 4570 + 60 K, log(g) = 1.7 + 0.1, [Fe/H] = -0.9 x 0.1, 
[a/Fe] = 0.7 + 0.1, and Vmicro = 0.6 + 0.1kms~!. We adopt the 
parameters from iSpec as our standard. Figure 2 compares a model 
spectrum generated using the iSpec parameters to the LBT/PEPSI 
spectrum. The model spectrum is a reasonable fit to the PEPSI data. 
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Table 1. Multi-band photometry measurements used in the construction of the SED for V723 Mon. Luminosities in each band are calculated assuming a nominal 


distance of d ~ 460 pc. 


Filter Magnitude c Fa [ergs s™! cm? A! ] ALalLo] Reference 

Swift UVM2 14.11 — 0.07 1.0 x 10714 0.15 This work 
SkyMapper u 10.30 0.01 2.5 x 107! 5.7 Onken et al. (2019) 
SkyMapper v 9.74 0.01 7.1 x 107! 18.1 Onken et al. (2019) 
Johnson B 924 0.05 1.3x 107? 36.5 Henden et al. (2018) 
Johnson V 830 0.04 1.7 x 1072 61.8 Henden et al. (2018) 
Sloan g’ 873 0.01 1.6x 10-2 49.1 Henden et al. (2018) 
Sloan 7.87 0.12 2.0 x 107!2 81.2 Henden et al. (2018) 
Sloan i? 749 0.14 1.9 x 1072 94.9 Henden et al. (2018) 
Sloan z' 7.36 0.05 1.1x 1072 71.8 Henden et al. (2018) 
2MASS J 6.26 0.03 9.8 x 107? 80.2  Skrutskie et al. (2006) 
2MASSH 5.58 0.03 6.7 x 107? 73.3  Skrutskie et al. (2006) 
2MASS K; 5.36 0.02 3.1 x 107? 43.7  Skrutskie et al. (2006) 
WISE W1 528 0.16 6.3 x 10714 14.0 Wright et al. (2010) 
WISE W2 5.10 0.07 2.2 x 10714 6.7 Wright et al. (2010) 
WISE W3 501 0.02 6.5 x 10716 0.5 Wright et al. (2010) 
WISE W4 4.82 0.04 6.0 x 1077 0.1 Wright et al. (2010) 


The spectroscopic parameters derived for the giant are summarized 
in Table 2. These estimates of the spectroscopic parameters do not 
consider the effects of veiling on the observed spectrum (see §4.4). 
Veiling introduces systematic uncertainties on these parameters and 
will lower the temperature estimate for the giant. 

The spectroscopic surface gravity is consistent with that inferred 
from the PHOEBE model in $4.3 (log(g)pHorBE = 1.63+0.06). Given 
the spectroscopic log(g) and the radius of the giant from §4.3, the 
spectroscopic mass is Mgiant,spec = 1.1 0.5 Mo, consistent with the 
mass of the giant derived from the PHOEBE model. The spectroscopic 
temperature is also consistent with that obtained from the SED fits. 
Based on the van Belle et al. (1999) temperature scale for giants, our 
temperature estimate is more consistent with a KO/K1 giant than the 
archival classification of GO (~6000 K) from Houk & Swift (2000). 
The absolute V-band magnitude (My =~ —0.3+0.1) is consistent with 
a luminosity class of III (Straizys & Kuriliene 1981). From single- 
star evolution, the spectroscopic measurement of log(g) suggests that 
the giant is currently evolving along the upper red giant branch. The 
giant has a luminosity larger than red clump stars, suggesting that it 
has not yet undergone a helium flash. 

We used the spectroscopic parameters in Table 2 and the luminos- 
ity constraint from the SED fit as priors to infer the physical properties 
of the giant using MESA Isochrones and Stellar Tracks (MIST; Dot- 
ter 2016; Choi et al. 2016). We used the isochrones package for the 
fitting (Morton 2015). We find that Mgiant.Mist = 1.07 + 0.24 Mo 
and Rgiant,MIST = 24.8 + 4.7 Ro. The age of the giant from the 


MIST models is MS. Gyr. These results are consistent with the 
properties of the giant derived from the spectra and the SED. 


4.2 Keplerian orbit models 


We fit Keplerian models both independently and jointly to the S12 
and G14 radial velocities using the Monte Carlo sampler TheJoker 
(Price-Whelan et al. 2017). The results for the four fits are summa- 
rized in Table 3. We first fit each data set independently as an elliptical 
orbit to verify that we obtain results consistent with the published 
results. We then fit the joint data set using either a circular or an ellip- 
tical orbit. In the joint fits we include an additional parameter to allow 
for any velocity zero point offsets between the S12 and G14 data. 
For the circular orbit we also set the argument of periastron w = 0. 
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Table 2. Properties of the red giant in V723 Mon (not accounting for veiling, 
see $4.4) 


Parameter FASMA iSpec 
Tar (K) | 4480450 4570460 
log(g) | 17202 1.7 £ 0.1 
[Fe/H] -1.1+0.1 -0.9+40.1 
[a/Fe] — 0.7 x 0.1 
Vmic (km s^!) 1.64 (fixed) — 0.64 € 0.07 
Vmac (km s^!) — 5.1 (fixed) 5.7 (fixed) 
vrot Sin i. (km s~!)* 19.1 x 0.4 
R (Ro) 24.9 + 0.7 
L (Lo) 173 +8 
M (Mo) 1.00 + 0.07 


* wot sin i varies with orbital phase 


Since the orbit is nearly circular even for the elliptical models, we 
derived the epoch of maximal radial velocity Tgy max instead of the 
epoch of periastron. We define phases so that Try max (BJD/TDB) 
corresponds to ¢ = 0.75, the companion eclipses at ¢ = 0.5 and the 
giant would be eclipsed at ¢ = 0. After doing a first fit for the ellip- 
tical models, we did a further fit with P, K, e, y and Try max fixed 
to their posterior values and further optimized w using least-squares 
minimization. 

The results of the fits are summarized in Table 3 and shown in 
Fig. 3. The fits to the individual data sets agree well with the published 
results, and the mass functions are well-constrained and mutually 
consistent. The elliptical models all yield a small, non-zero ellipticity, 
consistent with G14's arguments. We do find a small velocity offset 
of AV = 0.58 + 0.14 km s^! for the elliptical model and AV = 
0.84 + 0.21 km s^! for the circular model between the $12 and G14 
data. While the velocity residuals of the fits are small compared to 
K (x 1.7 km s^! versus 65 km s7!), they are large compared to 
the measurement uncertainties. Thus, while the RV curve is clearly 
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Figure 2. The LBT/PEPSI spectrum from 5400 À to 5600 A (black). A model spectrum using the atmospheric parameters derived from iSpec (Table 2) is 


shown in blue. 


completely dominated by the orbital model, Fig. 3 also shows that 
there are significant velocity residuals for both joint fits. The circular 
fit is dominated by a residual of period Po,,/2 and the elliptical fit is 
dominated by a residual of period P,,,/3. Fitting an elliptical orbit 
with a circular orbit will show a dominant Po,,/2 residual, but an 
elliptical orbit should not show a large Po;p/3 residual. This Pop /3 
residual to the fit of an elliptical orbit is likely the origin of the 
S12 hypothesis that the companion is an SB1 with this period. The 
residuals do not however, resemble those of a Keplerian orbit. We 
discuss this hypothesis and the velocity residuals in Appendix A. 


Binaries with evolved components (log(g) < 2.5) and orbital 
periods shorter than ~100 days are expected to have gone through 
tidal circularization and have circular orbits (e.g., Verbunt & Phinney 
1995; Price-Whelan & Goodman 2018). However, in the joint fits, 
the models with ellipticity are a better fit and have smaller RMS 
residuals than the circular models. Griffin (2014) carried out an F- 
test and noted that the ellipticity in their best-fit orbit for V723 Mon 
was significant, and concluded that it was very likely real. While we 
use the circular orbit for the PHOEBE models in §4.3, the differences 
in the mass function and semi-major axis compared to the elliptical 
orbit are small and have no effect on our conclusions. 


Independent of this issue with the origin of the RV residuals, the 
companion to the red giant in V723 Mon must be very massive given 
amass function of 1.72—1.74 Mo. The mass function itself is greater 
than the Chandrasekhar mass of ~1.4Mo, immediately ruling out a 
white dwarf companion. For an edge on orbit (see $4.3) and Mgiant = 
1 Meo, the companion mass is Mcomp ~ 3Mo and the semi-major 
axis is a ~ 100Ro. The Roche limits are approximately RL siant = 
0.29a = 30Ro for the giant and Rr, comp = 0.48a = 49Ro for the 
companion. Based on the radius estimate from the SED fits (§4.1), 
the giant is comfortably inside its Roche lobe (Rgiant/RL,giant = 
0.66) but in the regime where we should be seeing strong ellipsoidal 
variability due to the tidal deformation of the giant by the gravity of 
the companion (e.g., Morris 1985). Any stellar companion also has to 


be well within its Roche lobe or it would dominate the SED because 
its Roche lobe is significantly larger. Hence, we can be confident that 
we have a detached binary whose light is dominated by a giant that 
should show ellipsoidal variability and might show eclipses. 


4.3 PHOEBE and ELC binary models 


While it has been previously claimed that V723 Mon is a 
contact/semi-detached binary of the 8-Lyrae type, we find that this 
is very unlikely. As we just argued, both the giant and any compan- 
ion must lie well within their Roche lobes given the properties of 
the giant and the orbit. Additionally, the morphology of the light 
curve is inconsistent with those of detached and most semi-detached 
eclipsing binaries. Here we interpret the light curves as ellipsoidal 
variability and deduce limits on any eclipses of the companion for 
use in $4.5. 

We fit the ASAS V-band, KELT Rg -band and TESS T-band light 
curves (Figure 5) using PHOEBE 2.3 (Prša et al. 2016; Horvat et al. 
2018; Conroy et al. 2020). Since the companion appears to be dark 
and producing no eclipses, we fix it to be a small (R = 3 x 1076 Ro), 
cold (Te = 300 K) black body, use the simplest and fastest eclipse 
model (eclipse method-only. horizon) and do not include the 
effects of irradiation and reflection. We adopt the joint-circular RV 
solution from Table 3 (period, e = 0, systemic velocity, semi-major 
axis, and argument of periastron w = 0). We set the bolometric grav- 
ity darkening coefficient as B = 0.54, comparable to those obtained 
by Claret & Bloemen (2011) in the R-band for stars with T.g —4600 K, 
-] € [Fe/H] x -0.5 and 0 x log(g) < 2. We did not find signifi- 
cant differences in the final parameters when we varied 8 by +20% 
in the PHOEBE models. A 20% change in £ roughly corresponds to 
AT.g = 500 K or Alog(g) = 1.5. 

We initially performed trial fits to the KELT and TESS light curves 
simultaneously using the Nelder-Mead simplex optimization routine 
(Lagarias et al. 1998). The parameters from the trial fit were then used 
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Figure 3. The observed radial velocities for V723 Mon as a function of orbital phase, defined with the epoch of maximal RV at $ = 0.75 (top). The RVs from 
Strassmeier et al. (2012) [S12] are shown as black circles and the RVs from Griffin (2014) [G14] are shown as blue squares. The joint, circular (elliptical) RV 
fit from Table 3 is shown in red (purple). The two models closely overlap and are hard to distinguish. The residuals from both the circular and elliptical fit to the 
combined data are shown in the lower panels. The RV residuals are most likely a result of the ellipsoidal variability rather than the result of a triple system (see 


appendix §A.) 


to initialize a MCMC sampler with 20 walkers, which was then run for 
3000 iterations using the emcee (Foreman-Mackey et al. 2013) solver 
in PHOEBE 2.3. We fit for the binary mass ratio q = Meiant/ Mcomp, 
orbital inclination (7), the radius (Rgiant) and the effective temperature 
of the giant (Teff giant). We marginalize over the semi-major axis of 
the red giant (agian sini), the passband luminosities and nuisance 
parameters that account for underestimated errors in the KELT Rx- 
band and TESS T-band light curves. We adopt Gaussian priors on 
the effective temperature (Teff, giant = 4440 + 90 K) and the radius 
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(Rgiant = 22.2 + 0.8 Ro) based on the SED fits in $4.1. We adopt 
uniform priors of [70?, 90°] and [0, 1] for the orbital inclination and 
mass ratio respectively. 


The corner plot of the posterior samples for q, i, R giant and Terr, giant 
are shown in Figure 4 and the results of the best-fitting PHOEBE model 
are listed in Table 4. The errors in the parameters were derived from 
the MCMC chains. Our model is a good fit to the KELT Rx -band and 
TESS T-band light curves, shown in Figure 5. The PHOEBE model 
indicates that the orbital inclination of V723 Mon is nearly edge 


Table 3. Orbital Elements for V723 Mon 
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Parameter S12 G14 S12 + G14 S12 + G14 
Pow(d) 59.9358 + 0.0017 59.9391 + 0.0020 59.9394 + 0.0014 59.9398 + 0.0020 

K (kms™!) 65.483 + 0.068 65.209 + 0.117 65.360 + 0.081 65.150 + 0.125 
e 0.0152 + 0.0010 0.0186 + 0.0014 0.0158 + 0.0012 0 (fixed) 

w (°) 89.1 + 4.1 88.2 + 5.8 88.4 + 4.6 0 (fixed) 

y (km s7!) 1.81 + 0.053 3.02 + 0.07 1.95 x 0.07 1.88 + 0.10 


Coiant sini (Ro) 


71.567 x 0.081 
4096.196 x 0.696 


71.242 x 0.139 
4098.946 x 0.954 


77.425 + 0.096 
4096.529 + 0.773 


77.187 + 0.148 
4096. 168 + 0.038 


Trv max (BJD-2450000) 


RMS Residual (km s7!) 0.519 


0.969 0.876 1.17 


f (M) (Mo) 1.743 x 0.005 


1.721 x 0.009 


1.734 x 0.006 1.717 x 0.010 


Table 4. PHOEBE parameter estimates for the dark companion primary (DC) 
and red giant secondary (RG) 


Parameter DC RG 
Pos (d) 59.9398 (fixed) 
c (°) 0 (fixed) 
e 0 (fixed) 
y (km s!) 1.88 (fixed) 
asini (Ro)  25.8976'122! —— 77.178019 
i (?) 87.0717 
Tar (K) 300 (fixed) 4590 + 10 
R (Ro) 3x 1076 (fixed) 24.9 + 0.7 
4 0.328893 01076 
M (Mo) 3.04 + 0.06 1.00 + 0.07 


on I7 p ). The semi-major axis of the binary is agp sini = 
giant SIn i + ac sini = 102.7 + 1.3 Ro. The radius derived from the 
PHOEBE model (24.9 + 0.7 Ro) agrees well with those obtained from 
the SED fits and the MIST evolutionary models in $4.1. 

To verify the results of the PHOEBE models, we independently fit 
the KELT Rx-band light curve and the S12 radial velocities us- 
ing the differential evolution Markov Chain optimizer (Ter Braak 
2006) and the ELC binary modelling code (Orosz & Hauschildt 
2000). The optimizer was initialized using the parameters from the 
Nelder-Mead simplex optimization from PHOEBE and the optimizer 
was run for 10,000 iterations using 18 chains. We fit for the mass, 
radius and temperature of the red giant (Mgiant> Rgiant, Tet giant). 
the orbital inclination (7), the binary mass ratio (q), and the ob- 
served radial velocity semi-amplitude (K) assuming uniform priors. 
We included priors for the temperature (Teff giant = 4440 + 90 K), 
radius (Rgiant = 22.2 + 0.8 Ro), surface gravity(log(g) = 1.7 + 
0.2) and vyo¢ sing = 19.1 + 1 km s-l. From the ELC posteri- 
ors, we obtain Mgiant = 0.99+0:04 Mo, Tost giant = 457000 K, 


—0.05 
Rite = 24597 Ro, i = BESUDS q = 033031105. aud 


5 -1.6°’ 


K = 65.213 + 0.025 km s~!. This implies that the mass of the com- 
panion is Mcomp = XU Mo. The parameters obtained from 
ELC agree well with those obtained from PHOEBE. 


With the well-constrained Keplerian model for the RVs and the 


Reiant = 24.9707 Ro 


q=0.33+2:02 


Tese = 4590735 K 


Fan i q Test 


Figure 4. The corner plot for the posterior samples for q, i, Rgiant and Teff, giant 
in the best-fitting PHOEBE model to the KELT and TESS light curves. 


PHOEBE model for the ellipsoidal variations, we are able to directly 
determine the masses of the two components. The mass of the com- 
panion is 


(4) 


FM) say 22 Mo | f(M) IE 
Mcomp = = 


sin? i sin? i (1.7 Mo 1.3 

and from the PHOEBE models, we find that the red giant has a mass 
Mgiant =~ 1.00 + 0.07 Mo and the companion has a mass Mcomp ~ 
3.04 + 0.06 Mo. The reported errors are purely statistical and do 
not consider systematic effects (veiling, etc.) in the derivation of 
the binary solution. However, our results in §4.4 indicate that the 
veiling in the KELT and TESS light curves is minimal and should 
not significantly bias the mass estimates. In the absence of evidence 
for a stellar companion in the SED, spectra, or in eclipse, perhaps 
the simplest explanation for the mass of the companion (Mcomp ~ 
3.04 + 0.06 Mo) is that it is a non-interacting black hole in the 
“mass-gap” between 3 — 5 Mo (Özel et al. 2010; Farr et al. 2011). 


We discuss other scenarios in §5. The estimated mass of the red 
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Figure 5. The normalized ROAD B-band, ASAS, KELT and TESS light curves for V723 Mon as a function of orbital phase (defined with the epoch of maximal 
RV at $ = 0.75). The best-fitting binary models from PHOEBE are also shown. Second light from the veiling component (see $4.4) is considered in the models 
for the B (fgecond,B 6390) and V -band (Fecona,v~24%) data. 
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giant (Mgiant ~ 1.00 + 0.07 Mo) places it towards the lower end of 
measured red giant masses in the APOKASC catalog (Pinsonneault 
et al. 2018). Of the APOKASC sources with measured asteroseismic 
masses, only ~13% had masses lower than 1.0 Mo. The mass derived 
from the binary modelling is also very similar to the MIST estimate 
in $4.1. 

Given the configuration of the binary and the large size of the giant 
(^25 Ro), a companion with a radius of 1 Ro should be eclipsed 
for inclination angles i z 76°. Figure 6 shows the eclipses predicted 
for the ASAS, KELT and TESS light curves for main sequence stars 
with masses of 1 Mo, 1.5 Mo and 2 Mo. Any such stars would 
have produced relatively easy to detect eclipses. We will focus on the 
KELT limits because these data have higher S/N compared to ASAS 
photometry, and TESS only observed the eclipse of the red giant 
and not of the companion. At the expected phase of the eclipse, the 
RMS of the KELT data relative to the eclipse-free ellipsoidal (ELL) 
model is only 0.84% for phases from 0.46 to 0.54. If we bin the data 
0.01 in phase, the RMS of the binned data is only 0.38%. We will 
adopt a very conservative KELT band eclipse limit of 1%. Similarly, 
based on the Swift, ASAS and ROAD light curves, we will adopt 
eclipse limits of 10%, 3% and 2% for the UV M2, B and V bands, 
respectively. The Swift UV M2 light curve is shown in Figure 7. The 
UV M2 light curve appears somewhat variable, with one larger outlier 
on UT 2020-11-13 (¢ = 0.35). However, we do not see any evidence 
for an eclipse at the (conservative) level of ~10% at the expected 
phases 0.46 x $ x 0.54. While the TESS light curve does not cover 
the eclipse of the companion at ¢ = 0.5, it does cover the eclipse of 
the giant at ¢ = 0. For this eclipse, we estimate a conservative limit 
on the TESS band eclipse depth of ~0.3%. 


4.4 Veiling 


After the KELT Rx-band and TESS T-band light curves were jointly 
fit, we compared the resulting model V-band light curve to the ASAS 
observations, and found that the observed V-band variability ampli- 
tude was smaller than predicted by our model. Lower amplitudes 
at (usually) bluer wavelengths are generally attributed to an addi- 
tional source of diluting ("second") light. For X-ray binaries, it is 
also known as “disk veiling” due to additional light from accretion 
(e.g., Hynes et al. 2005; Wu et al. 2015). The amount of additional 
flux is usually characterized either by the ratio r of the veiling flux 
to the stellar flux or the veiling flux as a fraction of the total flux 
f =r/( +r). Interacting X-ray binaries are frequently observed to 
have veiling factors upwards of fgisk~50%, corresponding to r~1 
(see, for e.g., Wu et al. 2015). 

We fixed the system parameters to those from the joint fit to the 
KELT+TESS light curves and fit for additional flux (second light) in 
the V-band. The fit to the ASAS V-band light curve is consistent with 
an additional V-band flux (F2) amounting to fgecond,v = F2/(F2 + 
Feiant) = 24+2%. Figure 5 shows the PHOEBE model fit for the V-band 
after accounting for this additional flux. The redder ROAD Rc and 
Ie light curves agree well with the PHOEBE model without additional 
flux. The ROAD observations also agree with the PHOEBE model 
for the V-band with the additional flux (foecona,v = 24 + 2%) added. 
The discrepancy between the B-band light curve and a model without 
extra flux is even larger than for the V-band. Fitting the B-band data 
requires a larger second light contribution of foecond,B = 63 + 2%. 
The color of the extra light B — V = —0.2 + 0.1 is very blue even 
before applying any extinction corrections. This corresponds to stars 
with Tag ranging from ~13, 000 K to ~30, 000 K (Papaj et al. 1993). 
However, no hot star could have a rapidly rising SED from V to 
B and then drop rapidly from B to the UV M2-band. Stellar SEDs 
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are intrinsically broader in wavelength, like the model in Figure 1. 
Additionally, stars with these temperatures are easily ruled out by 
constraints from the SED and eclipses (see 84.5). 

Veiling also affects the spectrum of the star (“line veiling”) by 
reducing the observed depths of stellar absorption features. This is 
a well-known problem in interacting X-ray binaries (see, for e.g., 
Casares et al. 1993) and T-Tauri systems (Gahm et al. 2008). Veiling 
in spectroscopic observations is characterized through the fractional 
veiling 

V() 
r(4) = sa’ (5) 
where V (4) is the veiling spectrum and S(A) is the uncontaminated 
spectrum of the star (Casares et al. 1993). The ratio r (4) is related to 
the second light (F2) as 


F= r(A) Fgiant: (6) 


where Fgiant is the flux from the giant. The observed spectrum of the 
star is a function of the veiling factor Fi (4) = 1 4 r(4) 


r(4) = A(A)S(A) Frei (A), (7) 


where A(A) is a normalization factor. The monochromatic veiling 
factor 


EW,; 


—— m (8) 
EWobs 


Fyeil = 
is the ratio of the observed equivalent widths to those predicted from 
the standard/synthetic spectrum. 

We calculated the monochromatic veiling factors for various Ca 1 
and Fe1 absorption lines in the SES spectra and the red side of the 
PEPSI spectrum (A > 500 nm). For the Ca lines, we assumed an 
abundance of [Ca/H] = 0.02 + 0.02 derived from the PEPSI spec- 
trum. We generated a synthetic iSpec/SPECTRUM (Gray & Corbally 
1994) spectrum for the red giant using the atmospheric parameters 
in Table 2. Since the R ~ 220, 000 PEPSI spectrum has lower S/N at 
blue wavelengths, we use the mean veiling factors in the SES spectra 
for the absorption lines blue-ward of 500 nm. We use the standard 
deviations of the estimates from the individual SES spectra to es- 
timate the uncertainties. For the redder lines at > 500 nm, we use 
the veiling factors derived from the PEPSI spectrum at $ ~ 0.63 and 
assign errors of +0.10 in r(A). This method can have large systematic 
uncertainties (see Casares et al. 1993), but it provides a way to inde- 
pendently test the estimates from the dependence of the ellipsoidal 
variability amplitudes on wavelength. 

Figure 8 shows that the fractional veiling and second light as a 
function of wavelength r(A) rises steeply towards bluer wavelengths, 
with a power-law (4%) index of a = —5.0 + 1.5, steeper than that 
expected from the Rayleigh-Jeans law (a = —4). These results also 
indicate minimal contamination in the KELT Rx and TESS T filters 
that are used for our primary fits to the ellipsoidal variations. As a 
precaution, we ran PHOEBE models adding 10 — 20% extra flux in the 
KELT Rx band, and found that the masses of the red giant and dark 
companion are comparable to those obtained without considering 
veiling given the error estimates. The fractional second light in the 
V-band from the PHOEBE fits is comparable to that derived from the 
line veiling estimates. The PHOEBE estimate of the second light in the 
B-band is somewhat larger than that seen in the line veiling, however 
both estimates indicate a large contribution to the B-band flux. Based 
on the lack of eclipses in the light curves (see $4.5), we can also infer 
that this emission component has to be diffuse. 

Since the second light has a non-negligible contribution to the 
total SED of V723 Mon for 2 x 600 nm, the inferred properties 
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Figure 6. The observed ASAS, KELT and TESS light curves compared to various eclipsing models with main sequence companions to the red giant. The 
ellipsoidal model is shown in red. The depths of the eclipses at œ = 0.5 compared to the non-eclipsing ellipsoidal model is also shown. Eclipse constraints rule 
out single and binary stellar companions with Mgingle > 0.80 Mo and Mpinary > 1.51 Mo (see $4.5). 
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Figure 7. The Swift UVM2 light curve. The expected eclipse duration (~ 5 days) of the companion is highlighted in red. 


of the star from §4.1 can also be affected. We refit the SED of 
the giant after removing the estimated veiling component to obtain 
Teff giant ~ 4150+ 100K, Lgiant = 1539 Lo, Rgiant ~ 24.040.9 Ro 
and E(B — V) = 0.069 + 0.044 (Figure 9). The veiling component 
contributes ~12% of the total SED flux, with most of the flux in the 
bluer wavelengths. We perform iSpec fits to the PEPSI spectrum 
after truncating it to the redder wavelengths at A > 600 nm that 
are minimally affected by veiling. We obtain T.g = 4350 + 50 K, 
log(g) = 2.0+ 0.1, [Fe/H] = -0.7 + 0.1, [@/Fe] = 0.3 + 0.1, and 
Vmicro = 0.8+0.1 km s7}, From the parameters listed in Table 2, these 
differ by AT ot giant = 220 K, Alog(g) = 0.3, and A[Fe/H] = 0.3. In 
contrast, the parameters derived from the PHOEBE fits (Table 4) do 
not change significantly if we change Ter giant- Because the models 
of the veiled SED are driven to larger, cooler stars, models forcing 
the star to be more compact are even more disfavored than in the 
unveiled models. The best model of the veiled SED has x? = 14.2, 
rising to 33.9 for Rgjant = 20 Ro and to 64.0 for Rgiant = 18 Ro. 


The origin of this veiling component remains unclear although it 
is clearly non-stellar in nature. However, the morphology of the veil- 
ing component is broadly compatible with the spectra of advection 
dominated accretion flows (ADAF). ADAF spectra can be described 
by contributions due to synchrotron emission, Compton scattering 
and Bremsstrahlung radiation (Quataert & Narayan 1999). The most 
luminous feature in the ADAF models come from the synchrotron 
peak and for stellar mass black holes, this peak falls in the optical 
wavelengths (Quataert & Narayan 1999). In their ADAF models for 
quiescent black hole binaries (Mg = 6 Mo) with low accretion rates 
(log(M/Meaa) ^ —4), Esin et al. (1997) show that the SED peaks 
at optical wavelengths and rapidly decays at both longer and shorter 
wavelengths. For V723 Mon, detailed ADAF models are necessary 
to determine whether the veiling component is related to accretion 
flows around the dark companion, but such analysis is beyond the 
scope of this paper. 


4.5 Limits on Luminous Stellar Companions 


We can constrain the presence of luminous companions using either 
the SED or the absence of eclipses. The limits using only the SED 
will be weaker than those using the eclipses, but are also independent 
of any knowledge of the inclination. For the SED constraints, we 
require that the companion contributes less than 100%, 60% and 
20% of the light in the UV M2, B and V bands respectively. The B 
and V band limits correspond to the estimated veiling source from 
$4.4 and are conservative since the SED of the veiling light appears 
to be inconsistent with a star, and the UV M2 band limit is simply the 
total observed flux because we lack any constraint on the amount of 
veiling for this band. For the eclipse constraints, we require that the 
companion contributes less than 10%, 3%, 2% and 1% in the UV M2, 
B, V and R bands based on the eclipse models in $4.3. While we did 
not use the UV M2 band in the SED fits because the SED models 
cannot account for any possible chromospheric emission from the 
rapidly rotating giant, there is no issue with using them to obtain 
these limits. 

To provide models for the relationships between luminosity, tem- 
perature and mass, we sampled stars from [Fe/H] = —1 PARSEC 
(Bressan et al. 2012) isochrones with ages from 1 to 10 Gyr, and we 
logarithmically interpolated along each isochrone to sample more 
densely in mass. We considered models in which the companion is 
either a single star or two stars. In the latter case, we considered all 
combinations of two stars on the same isochrone. The age of the stars 
is the principal variable leading to changes in the mass limits, but 
the differences are not large. Two star models will generically allow 
higher limiting masses for two reasons. First, since luminosity is a 
steeply rising function of mass, dividing a single star into two lower 
mass stars of the same total mass leads to a lower total luminosity. 
Second, the two lower mass stars are also cooler on the main se- 
quence, so they produce less blue light per unit luminosity than the 
single star, which also allows the total mass be larger for the same 
constraints on the bluer fluxes. 

We start with the weaker but inclination independent limits from 
the SED, where Figure 10 shows the mass limits as a function of 
stellar age. For a single star, the maximum allowed mass is 1.29 Mo 
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Figure 8. The fractional veiling r (4) (Equation 5) [left] and second light (Equation 6) [right] calculated using Ca1 and Fer absorption lines as a function of 
wavelength. The fractional second light in the BV -bands as calculated from PHOEBE are shown as filled diamonds. Power-law fits are shown as blue dashed 
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Figure 9. The best-fitting SED model for V723 Mon after correcting for the 
veiling component (red triangles). The uncorrected SED from Figure 1 is 
also shown (dashed line). Note that the SED of the veiling component is very 
different from the SEDs of stars, as can be seen by comparing to Figure 1. 


(log(age) = 9.49, L = 15.5 Lo, T = 5430 K), and it occurs where 
the star is just starting to evolve off the main sequence because the 
initial drop in temperature weakens the constraints more than the 
rise in luminosity (see Figure 1). There are two maximum masses for 
models with two stars. The lower age peak corresponds to pairing a 
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slightly evolved star with a lower mass main sequence star. The total 
mass of 1.98 Mo comes from combining a 1.25 Mọ (L = 12.0 Lo, 
T = 5510 K) star with a 0.77Mo (L = 0.5 Lo, T = 5900 K) star at 
log(age) = 9.45. The higher age peak comes from combining two 
stars of the same mass. The total mass is again 1.98 Mo and the two 
components each have M = 0.99 Mo (log(age) = 9.80, L = 8.5 Lo, 
T = 5430 K). Fig. 1 shows the SEDs of these three models as 
compared to that of the giant, and they are nearly identical. 


For the BH candidates LB1 and HR 6819, the stellar companions 
seemed "dark" because they were hot. This cannot be the case here 
because of the tight UVM2 luminosity limit of ALa « 2 Lo (after 
correcting for extinction). Models normalized by this luminosity only 
have total luminosities of 4 Lo, 5 Le and 75 Lo, for temperatures of 
10^ K, 3x 10^ K and 10? K respectively, much too low for a massive 
He star companion. For example, a Helium star with T ~ 10? K will 
have M ~ 8 Mo and L ~ 90,000 Lo (Gráfener et al. 2012). For any 
hot companion with a luminosity approaching that of the giant, we 
would also expect to see the effects of irradiation of the giant on its 
light curve and phase dependent spectra, but no such perturbations 
are seen. 


Given the size of the orbit and the size of the giant, compan- 
ions with radii similar to those allowed by the SED will be eclipsed 
provided the inclination angle isi = 76°. While the systematic uncer- 
tainties on the estimated inclination of ERU may be moderately 
larger, the light curve shapes at inclinations anywhere approaching 
this upper limit for seeing eclipses are grossly incompatible with the 
observations. Fig. 10 also shows the limits on masses using the flux 
limits on the companion required to avoid visible eclipses, and they 
are stronger than those from the SED as expected. The biggest change 
is that they eliminate the “bumps” associated with stars starting to 
move off the main sequence because they more directly constrain 
the size of any companion. The improvements are otherwise modest 
because the (blue) luminosities are such strong functions of mass 
that moderate changes in the limits on the luminosity produce very 
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Figure 10. Limits on the mass of a companion that is comprised of either one 
(black) or two (red) luminous stars based on either the lack of KELT eclipses 
(solid) or the SED (dashed) as a function of the age of the stars. From these 
limits, the maximum allowed mass for a single star is 1.29 Mo (SED) and 
0.80 Mo (eclipses). For two stars, the maximum allowed mass is 1.98 Mo 
(SED) and 1.51 Mo (eclipses). 


modest changes in the limits on the mass. The single star limit is now 
0.80Mo6 and the two star limit is 1.51 Mo where the two star limit 
is always weakest for the stars being twins. These limits are smaller 
than the binary mass function for this system (see $4.2). 


4.6 Balmer Ha and Hf emission 


We find that the Balmer Ha and Hf lines appear to significantly vary 
with phase (Figure 11) and this is unusual for a red giant. Given our 
phase convention, the dark companion will be eclipsed by the giant 
at @ = 0.50 with an eclipse duration of tec] = 5 days. If a companion 
is responsible for Balmer emission, the subtraction of a template 
spectrum near phase ¢@ = 0.50 should isolate its contribution. To 
explore this, we subtracted the SES spectrum of the red giant near 
$ ~ 0.5 from the remaining spectra with S/N > 30. 

Figure 11 shows the changes in the Ha, Hf, Car46439 and 
Ca146463 lines with orbital phase. The left panel shows the con- 
tinuum normalized spectra, the middle panel shows the equivalent 
widths and the right panel shows the residuals relative to the spec- 
trum observed when the companion should be in eclipse. A first point 
to note is that all four lines show a phase-dependent change in equiv- 
alent width which mirrors the ellipsoidal variability - this effect is 
well-known (e.g., Neilsen et al. 2008) and further confirms the origin 
of the variability ($4.3). The Balmer absorption lines appear to be 
blue shifted by ~12 km s7! (Figure 12). This is also seen in the spec- 
tra close to conjunction at @ = 0.5 (Figure 11). However, we do not 
see a similar shift in the other photospheric lines. The second point 
to note is that while photospheric absorption lines like Car 46439 
and Ca146463 lines are cleanly subtracted, the Ha and Hf emission 
clearly varies with orbital phase. We do not see the Cau H and K 
lines in emission, so the changes in the Balmer lines are unlikely to 
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be caused by chromospheric activity. The Balmer emission could be 
caused by mass loss from the red giant through a stellar wind. 

The typical Gaussian FWHM of the Ha and H£ emission profiles 
is ~290kms~!. The median equivalent width of the residual Ha 
and Hf lines is EW(Ho) = 1.79 + 0.02 and EW(Hf) = 2.72 + 
0.02 respectively. We can convert the Ha equivalent width to the 
flux at the stellar surface using Fg, = Fe EW(Ha) (see for e.g., 
Soderblom et al. 1993; González Hernández & Casares 2010). Fe is 
the continuum flux which we derive using Hall (1996) as log( Fc) = 
7.538 — 1.081(B — V)o. For V723 Mon, we find (B — V)g = 0.85 + 
0.06 mag, using the APASS DRIO photometry and the E(B — V) 
from the SED fits. We obtain log(F.) = 6.61 + 0.06 and Fg, = 
7.3 + 1.0 ergs cm"? sl. Normalizing the Ha flux to the bolometric 
flux (Le. Rya = Fga/o0 T^), we obtain log(Rya) = —3.53 + 0.06. 
log( Rp) is usually compared to the Rossby number Ro = Prot/Te ~ 
2.84, where Te ~ 21.1 d based on Noyes et al. (1984) and assuming 
that the giant’s rotation is fully synchronized with the orbital period. 
At this value of Ro, V723 Mon has a value of log( Rp) higher than 
any of the chromospherically active single stars from López-Santiago 
et al. (2010) (see figure 5 in González Hernández & Casares 2010). 
This also indicates that the observed Ha emission is not just from 
chromospheric activity. 

Another argument against chromospheric activity is that the 
changes in the Balmer lines with phase are coherent over the 
—3.5 years spanned by the SES data. If the Balmer emission was 
dominated by chromospheric activity, we would expect the structure 
to change with time as the spot patterns evolve. At RV quadrature 
(@ 2 0.25 and 9. = 0.75), the emission profiles of the Balmer lines 
resemble P-Cygni profiles with both absorption and emission compo- 
nents. This is illustrated in Figure 13. In comparison, the Ca 1 46439 
line is cleanly subtracted and does not show a similar correlation with 
the phase of the binary orbit. P Cygni profiles tend to be associated 
with mass outflow and stellar winds. At RV minimum (¢ = 0.25), 
the absorption component in the Balmer lines is red shifted whereas 
at RV maximum (¢ « 0.75), the absorption is blue shifted. In the 
rest-frame of the giant, the peak of the Ha emission component is 
relatively stationary at —35 km s71. The median separation between 
the absorption and emission components in Ha is AV~115 km st, 
however this appears to vary with orbital phase. The separation is 
largest for 0 < $ < 0.25, with AV~148kms7!, and drops there- 
after. The smallest separation between the two components is seen 
at phases 0.5 < $ < 0.75, with AV~104kms~!. Much like with the 
SES spectra, we also see clear variability in the Ha line profiles in 
the HIRES spectra (Figure 14), and the asymmetry in the HIRES line 
profiles also reverses after $ = 0.5. 

Remarkably, near conjunction (¢ ~ 0.5) there is very little Balmer 
emission (Figure 13). Both the Balmer and photospheric lines are 
modulated with the ellipsoidal variations (Figure 11), however, near 
$ = 0.5, the equivalent width of the Balmer lines increases abruptly, 
signalling a dramatic drop in Balmer emission. A similar EW increase 
is not seen in the photospheric lines. This feature is coincident with 
the eclipse of the unseen companion by the red giant and has the 
expected duration. The MODS spectra were taken around the eclipse 
of the dark companion by the red giant at @ ~ 0.5 and the depths 
of the MODS Balmer Ha and Hf absorption features also deepen 
during the eclipse. 

While we have clear evidence to show the presence of Balmer 
emission that is correlated with the orbital motion of the putative 
black hole (84.6), the exact origin of the Balmer emission is unclear. 
One explanation for this additional luminosity and the variability in 
the Balmer lines is from an accretion flow. Alternative explanations 
for the Balmer emission can come from a stellar wind and/or mass 
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Figure 11. SES line profiles for the Balmer Ha, H8, Car46439 and Ca 146463 lines as a function of orbital phase (left). The middle column shows the 
equivalent width normalized to ø ~ 0 as a function of phase. A simple sinusoidal model (red) is fit to the variations in the Ca 146463 line to illustrate the 
variations in the equivalent width due to ellipsoidal variations. The shaded gray region illustrates the duration of the eclipse (~5 days). The right column shows 
the residuals of the line profiles after a template for the giant is subtracted. 
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Figure 12. LBT/PEPSI line profiles for the Balmer Ha, Hf, Ca146439 and Ca 146463 lines (black). A model spectrum using the atmospheric parameters 
derived in $4.1 is shown in red. The blue lines show the velocity offset of the Balmer absorption lines (^12 km s7!) from the rest frame of the giant. 


outflow at the inner Lagrangian point (L4). At ¢ = 0, the inner La- 
grangian point is directed toward the observer. This also coincides 
with the deeper minimum in the light curve because the surface 
gravity and brightness is smallest at Lı (Beech 1985). Conversely, 
at @ = 0.5, L4 points away from the observer. It is possible that 
the Balmer emission is associated with the photoionization of matter 
streaming through the inner Lagrangian point. The nearly-constant 
velocity offset of the emission peak from the secondary is also consis- 
tent with this interpretation. However, this scenario requires a source 
of photoionization. The neutron star binary 1FGL J1417.7—4402 also 
has a persistent Ha line with a complex morphology (Strader et al. 
2015; Swihart et al. 2018). Instead of an accretion disk, the authors 
attributed this behavior to the interaction between the magnetically 
driven wind from the secondary and the pulsar wind. Balmer photons 
originating from an interbinary shock or the wind from the secondary 
will have a velocity offset from the secondary. This is also seen in 
V723 Mon. However, the Ha line in 1FGL J1417.7—4402 is sig- 
nificantly different from what is observed in V723 Mon as it has a 
double-peaked emission profile that is observed directly in the spec- 
tra without the need for template subtraction. Swihart et al. (2018) 
also did not find independent evidence of a disk through their light 
curves, whereas we see clear evidence of second light in the B and 
V-band light curves for V723 Mon. 


In summary, we see evidence of broad Balmer emission in the 
spectra. The Balmer line profiles prior to when the unseen compan- 
ion is eclipsed by the giant (0 € ¢ x 0.46) resemble P-Cygni profiles, 


with a blue shifted emission component and a red shifted absorption 
component in the rest-frame of the giant. When the unseen com- 
panion is behind the giant at ø ~ 0.5, both the Balmer components 
disappear and the Balmer absorption features from the giant become 
deeper. When the unseen companion re-emerges after the eclipse at 
phases 0.54 < $ < 1, we see a clear change in the P-Cygni-like 
line profile, with both the absorption and emission components blue 
shifted. However, the absorption component is blue shifted more 
than the emission and the P-Cygni profile is reversed. We also see 
significant changes in the Balmer line depths at $ = 0.5 when the 
companion is eclipsed. The exact origin of the Balmer emission, be 
it from an accretion disk around the black hole, a stellar wind, mass 
flow into an inner Lagrange point in the binary, or a combination of 
these, remains unclear. 


4.7 X-ray detection 


In the individual Swift XRT exposures, we only detected a source 
above the background in the exposures taken on 2020-11-05 (9 = 
0.2), 2021-01-21 (9 = 0.5) and 2021-02-20 (9 = 0.0) with 4.2 + 2.2, 
6.8 + 2.9 and 3.7 + 2.2 background-subtracted counts in the 0.3- 
2.0 keV energy range, respectively. In the merged Swift XRT data, 
we detect a source with 18.5 + 4.9 background-subtracted counts in 
the 0.3-2.0 keV energy range, corresponding to an aperture cor- 
rected count rate of (1.4 + 0.4) x 107? counts/second. Assum- 
ing that E(B — V) = 0.086 from $4.1, and using the relation- 
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Figure 13. The right panels show the template subtracted Ha (red dot-dashed), HS (blue dashed) and Ca126439 (black) line profiles at various orbital 
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configurations illustrated in the left panels. The size of the compact object is simply chosen to make it easily visible. 
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Figure 14. The HIRES Ha line profiles as a function of orbital phase. The 
template SES spectrum at $ = 0.5 is shown in green. The variability in the 
line profiles is the same as that seen in the SES spectra (see Figure 11, $4.6). 


ship between reddening and atomic hydrogen column density from 
Liszt (2014), the Galactic column density towards V723 Mon is 
Nu,v723 = 7.1 x 1029 cm-2, This is considerably smaller than 
the total column density in the line of sight towards V723 Mon 
(NH,Los = 3.42 x 10?! cm2). Assuming an absorbed power-law 
with a photon index of 2 and the Galactic column density esti- 
mate Ny,v723, the Swift XRT count rate corresponds to an ab- 
sorbed flux of (3.0 + 0.8) x 10714 ergs cm"? s^! or an unabsorbed 
flux of (4.1 + 1.0) x 10-4 ergs cm 2 s^! in the 0.3-2.0 keV en- 
ergy range. At the Gaia EDR3 distance, we obtain absorbed and 
unabsorbed X-ray luminosities of (7.6 + 2.0) x 1029 ergs s^! and 
(1.0 + 0.3) x 109? ergs s~!, respectively. If we use a 0.3-10.0 keV 
energy range, we obtain a similar number of counts as there is very 
little emission >2.0 keV?. We also derive a limit on the hardness 
ratio! of —0.21 which indicates a relatively soft emission spectrum 
that is consistent with the lack of hard X-ray emission above 2.0 keV. 
After grouping the XRT observations by orbital phase, we calcu- 
late the out-of-eclipse and in-eclipse absorbed X-ray luminosities as 
(6.12:2.2) x 1029 ergs s! and (9.343.1)x 1029 ergs s^! respectively. 
The out-of-eclipse and in-eclipse X-ray luminosities are consistent 
given the measurement uncertainties. 

In the XMM-Newton observation, V723 Mon is not detected above 
the background, so we derive an upper-limit on Lx. We find a 3c 
upper-limit count rate of 4.17x10^? counts/second in the 0.5-8.0 
keV band within a 60" aperture (consistent with the off-axis PSF 


15 We derive a 30 upperlimit of 1.2 x 107^ counts/second to the 2.0-10.0 
keV count rate, which corresponds to an absorbed flux of 1.4 x 107!*ergs s^! 
and a luminosity of 3.6 x 1072 ergs s^! 

1^ Here the hardness ratio is defined as HR = (H — S)/(H + S) where H 
is the number of counts in the 2.0-10.0 keV energy range and S is the number 
of counts in the 0.3-2.0 keV energy range. 
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at the source position). Assuming an absorbed power-law with a 
photon index of 2 and the Galactic column density of Ng y723, we 
find an upper-limit absorbed flux of 1.0x10~!3 erg cm"? s7! and 
an unabsorbed flux of 1.1x107 P? erg cm? s^!. Adopting the Gaia 
EDR3 distance, the latter quantity yields an X-ray luminosity upper- 
limit of 2.7x1099 erg s7! in the 0.5-8.0 keV band. 

The Swift XRT estimate of the X-ray luminosity is comparable to 
the Lx ~ 1030 — 103! ergs s7! observed for quiescent X-ray binaries 
(Dincer et al. 2018) and the Lx ~ 1029 — 1030 ergs s^! observed 
for chromospherically active RS CVn systems (Demircan 1987). If 
the X-ray luminosity originates from an accretion disk, the putative 
black hole accretes at a very low luminosity of -107? Legg. There 
is observational evidence which indicates that the X-ray luminosity 
from quiescent black holes is significantly fainter than that from neu- 
tron stars (see for e.g., Asai et al. 1998; Menou et al. 1999). While 
debated, this has been attributed to advection-dominated accretion 
flows (ADAF; see, for e.g., Narayan & Yi 1995) and the existence of 
event horizons for black holes. However, given the rapid rotation of 
the giant, some of the observed X-ray luminosity may originate from 
the giant's chromosphere (Gondoin 2007). The X-ray spectra of X- 
ray binaries are generally harder (hotter), with average temperatures 
kT » 5 keV, compared to stellar coronae that have lower tempera- 
tures, with kT « 2 keV (Kong et al. 2002). Given the results of the 
Swift XRT analysis, it appears that the X-ray emission is relatively 
soft, which appears to be consistent with a chromospheric origin. 
However, to fully characterize this X-ray emission, followup X-ray 
spectra of significantly better S/N are necessary. 

We can estimate the mass loss rate from the giant as 


LR 
M 24x 107 P gg PELO x 8 x 107! Moyr!, (9) 


(Reimers 1975) where L, R and M are in solar units and nr = 
0.477 (McDonald & Zijlstra 2015). The velocity of the wind from 
the giant is assumed to be the escape velocity, which is Vwing % 
Vese ~ 120km s^l. For a scenario where the black hole accretes 
mass through the stellar wind, Thompson et al. (2019) approximated 
the amount of material gathered at the sphere of influence of the 
black hole as 


] 2 
: M GM 
Mace ~ n | 22E 
(4ra. ) 2 
giant wind (10) 


2 -4 
Mpu | Vwind ) 
3Mo] \ 120km s7! 


~7x107!! Mo yr! Wis 


where Mọ = M/ 107? Mo yr. For radiatively efficient accretion 
onto the black hole, the accretion luminosity can be approximated by 


Dace ~ 0.1 Mace c? = 83 Lo. D 


For radiatively inefficient accretion, the luminosity can be much lower 
(Narayan & Yi 1995). This estimate of the accretion luminosity is 
much larger than the observed X-ray luminosity of this system ($4.7), 
but closer to the luminosity of the veiling component (Lace ~% 4Lyeil). 


5 THE NATURE OF THE DARK COMPANION 


Given the observed properties of the system and the modeling results 
from §3 and $84, we next systematically discuss the nature of the 
companion. We first discuss the uncertainties in the mass of the 
companion and then systematically consider the possible single and 
binary possibilities for its composition. 
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Ultimately, our knowledge of the mass of the companion is de- 
termined by how well we can constrain the properties of the giant. 
Estimates of the mass through either the radius and gravity or stellar 
evolution models give relatively crude limits of 1.1 + 0.5 Mo and 
1.1 € 0.2 Mo, respectively (84.1). If the giant has a degenerate core, 
which is likely to be true because it is more luminous than the red 
clump, then there is a very firm lower limitto its mass. The luminosity 
of a giant with a degenerate core of mass Mcore is 


7 
L Moore 
= 67 , 12 
A " (ore V9 


(Boothroyd & Sackmann 1988), so we must have Maiant > Moore = 
0.33 Mo for L ~ 170 Lo. This implies a lower bound for the com- 
panion mass of Mcomp 2 2.26 Mo, above the mass of the most 
massive neutron star observed and larger than the limits for single 
and binary stellar companions in $4.5. 

The stronger limits on the mass come from the PHOEBE mod- 
els of the ellipsoidal variability in $4.3. Ignoring inclination, ellip- 
soidal variability constrains the mass ratio q because (to leading 
order) the amplitude of the Po,,/2 term depends on e; ~ (R/a)?/q 
while the amplitudes of the Pop and Porb/3 terms depend on 
(Rgiant/a)*/q (e.g., Morris 1985; Gomel et al. 2020). Since the semi- 
major axis is determined by the period, mass function and mass ratio, 
a? œ P? f(1+q)?, the ellipsoidal variability can constrain both Rgiant 
and q. This also means that the radius of the giant is an important 
independent constraint. Given the period, mass function and ampli- 
tude €», the mass of the giant Mgiant = fq(1- q)? œ RÀ. simply 


giant 
scales with the radius of the giant. The mass of the Conipatübhs 


Mcomp = f (1+ g’, is less dependent on the radius because the mass 
ratio q is small. We can verify this correlation using PHOEBE models 
with the radius of the giant fixed to Rgiant = 18, 20, 22, 24, and 26 Ro. 
The mass of the giant increases monotonically and fairly rapidly with 
radius, Mgiant = 0.42, 0.55, 0.72, 0.91, and 1.12 Mo, while the mass 
of the companion increases more slowly, Mcomp z 2.38, 2.55, 2.75, 
2.95, and 3.16 Mo, as expected. In fact, the numerical results al- 
most exactly track the expected scalings from holding e» fixed while 
varying the radius of the star. 

The PHOEBE models on their own found a radius of Rgiant = 24.9 + 
0.7 Ro, which agrees well with the independent results from the SED 
fits. Without the correction for veiling, these fits gave 22.2 + 0.8 Ro 
and with the correction for veiling they gave 24.0 + 0.9 Ro. These 
fits use a much more complete SED model than used for the Gaia 
DR2 estimate of ~20 Ro, but even for this smaller estimate the 
companion mass is Mcomp — 2.6 Mo. Since the radius estimates 
from the SED and the ellipsoidal variability agree, we will proceed 
under the assumption that the resulting mass estimates of Mgiant = 
1.00+0.07 Mo and Mcomp = 3.04+0.06 Mo are essentially correct. 

We next consider all the possible scenarios for the composition of 
the companion, with Table 5 providing a summary. If the companion 
is a single object, the options are a star, a white dwarf (WD), a neutron 
star (NS) or a black hole (BH). A star is ruled out by the eclipse and 
SED limits from §4.5, with Mcomp < 0.8Mo based on the lack of 
eclipses. A WD is ruled out because its mass would exceed the Chan- 
drasekhar mass limit of 1.4Mo. A NS is unlikely, since a companion 
mass of Mcomp ~ 3.04 Mo is still slightly larger than the maximum 
passable mass of a neutron star (M ~ 3 Mo; Lattimer & Prakash 
2001). Additionally, the mass is well above the maximum observed 
masses of M =~ 2.01 +0.04 Mo and M = 2. 14 a Mo found for the 
neutron stars PSR JO348+0432 and MSP J0740+6620, respectively 
(Antoniadis et al. 2013; Cromartie et al. 2020). The mass of the dark 
companion is slightly larger than the M ~ 2. 59+0.09 Mo mass-gap 


—0.08 
compact object in the LIGO/VIRGO gravitational wave merger event 
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GW190814 reo et al. 2020), and it is comparable in mass to the 
Mpy ~ 3. bee 7 Mo non-interacting black hole identified around the 
red giant MASS J0521565844359220 by Thompson et al. (2019). 
The mass estimates for V723 Mon, due to its ellipsoidal variability, 
are much tighter than for 2MASS J05215658+4359220. Thus, the 
simplest explanation of observed V723 Mon system properties is 
that the binary companion is a black hole near the lower end of the 
mass gap. 

There are many more possible scenarios if the companion is a 
binary. The orbit of such a binary has to be fairly compact (see 
Appendix A), but not too compact, or the system would merge too 
rapidly (see below). While a single phase of common envelope evo- 
lution can likely lead to the simple binary models, any of these triple 
solutions likely requires a more complex evolutionary pathway which 
we will not attempt to explore here. Based on our analyses in $4.5, 
the dark companion cannot be a stellar binary given that the two star 
mass limit is Mpinary < 1.5 Mo from the eclipse constraints (Figure 
10). The companion must contain at least one compact object. 

Combining a star with a compact object seems unlikely. With the 
stellar mass « 0.8 Mo based on the lack of eclipses, the compact 
object mass must be > 2.2 Mo. A WD is ruled out because this 
exceeds the Chandrasekhar limit. A NS is possible, but the mass is at 
or above the maximum observed NS masses. A BH is also possible, 
but a single more massive BH seems far more plausible than a BH 
with a NS-like mass combined with a star to form the inner binary 
of a triple system. 

A double WD binary requires that both WDs are very close to the 
Chandrasekhar limit or above it. Given that such massive WDs are 
quite rare (Tremblay et al. 2016), putting two of them into such a 
system is unlikely even if the masses can be kept just below the limit. 
Combining a WD and a NS is allowed, and the NS mass is in the 
observed range if the WD massis » 0.9Mo. Combining a WD with a 
BH is possible, but has the same plausibility problems as combining 
a star with a BH. A double NS binary is feasible. They would need 
to be of similar mass, since a mass ratio > 0.67 is needed to keep 
the lower mass NS above the theoretical minimum of 1.2 Mọ (Suwa 
et al. 2018). Combining a NS with a BH and a double BH binary 
seem implausible because they both require BH masses in the range 
observed for neutron stars. 

An additional consideration for any compact object binary model 
for the companion is its lifetime due to the emission of gravitational 
waves. For an equal mass, 3 Mo binary in a circular orbit with 
semi-major axis a;,, the merger time is 


3 Mo 

tmerge = 2.2 x 107 years (a EA (13) 
If the age of the system must be » 1 (10) Gyr to allow time for the 
red giant to evolve, then aj, > 2.6 Ro (4.6 Ro). Unfortunately, a 
binary with such a long life time is too weak a gravitational wave 
source to be detected by the Laser Interferometer Space Antenna 
(LISA, Robson et al. 2019). As discussed in Appendix A, dynamical 
stability requires aj; < 31Ro (significantly inside the Roche lobe 
radius of the companion, ~49 Rọ), so a long-lived, dynamically 
stable binary is possible for 4 Ro S aij S 31 Ro. This does not 
guarantee secular stability, and for much of this range of semi-major 
axes we would also expect to see dynamical perturbations of the outer 
orbit and additional tidal interactions (see Appendix A). 

In summary, while it is difficult to rule out more complex scenarios 
where the companion to the giant is a binary consisting of at least one 
compact object, the simplest explanation is that the dark companion 
is a low-mass black hole. This would make V723 Mon a unique 
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Table 5. Comparison of the possible scenarios involving stellar companions (*), white dwarfs (WD), neutron stars (NS) and black holes (BH) that can explain 
the nature of the dark companion. We assumed Mcomp = 3.04 + 0.06 Mo from 84.3 (Table 4). A */" indicates that the scenario is possible, a *?' indicates that 
while the scenario is technically possible, it is very unlikely and a ‘X’ indicates that the scenario is ruled out. The simplest explanation is that of a single low 


mass black hole, indicated with ‘V/V’. 


Dark Companion Possibility Comment 
Single Star x Ruled out by SED/eclipse limits from $4.5 (M. x 0.8 Mo). 
Single WD x WD will exceed Chandrasekhar limit (Mwp > 1.4 Mo). 

Single NS 2 Requires an extreme NS equation of state. 
Single BH Vv Simplest explanation. 
Star + Star x Ruled out by SED/eclipse limits from §4.5 (Mpinary S 1.5 Mo). 
Star + WD X For M. x 0.8 Mo (84.5), the WD mass exceeds Chandrasekhar limit. 
Star + NS y For M. « 0.8 Mo, the NS mass exceeds 2.2Mo. 
Star - BH ? BH mass is even lower than with no star. 
WD + WD x Both WD components near or above Chandrasekhar limit. 
NS + WD P4 NS mass is in the observed range if Mwp > 0.8 Mo 
BH + WD ? BH mass is even lower than with no WD. 
NS + BH x The BH must have a NS-like mass. 
NS + NS v Both NS components should have Mys z 1.2 Mo, so dinner Z 0.67. 
BH + BH x The BHs have NS masses. 


system, as it would contain both the lowest mass BH in a binary 
and be the closest known black hole yet discovered. The low X-ray 
luminosity (L/ Lega 107?) of this system likely also suggests a black 
hole companion, as quiescent black holes are known to be X-ray faint 
(see $4.7). 


6 CONCLUSIONS 


The nearby (d ~ 460 pc), bright (V ~ 8.3 mag), evolved red giant 
(Toff giant ~ 4440 K, Lgiant ~ 173 Lo) V723 Mon is in a high mass 
function, f(M) = 1.72+0.01 Mo, nearly circular binary (P = 59.9 d, 
e = 0) with a dark companion of mass Mcomp = 3.04 + 0.06 Mo. 
V723 Mon is a known variable that had been typically classified as an 
eclipsing binary, but the ASAS, KELT and TESS light curves indicate 
that is in fact a nearly edge-on ellipsoidal variable ($4.3, Figure 5). We 
do not see any continuum eclipses in these light curves ($4.5, Figure 
6). Using the binary mass function and sini ~ 1 ($4.3), itis clear that 
V723 Mon has a companion of mass 2.6 Mọ < Mcomp < 3.6 Mo 
for 0.5 Mo < Mgiant < 1.5 Mo. We modeled the light curves with 
PHOEBE using constraints on the period, radial velocities and stellar 
temperature to derive an inclination of i = BLU MS a mass ratio 
of q = 0.33 + 0.02, a companion mass of Mcomp = 3.04 + 0.06 Mo, 
a stellar radius of Rgjant = 24.9 + 0.7 Ro, and a giant mass of 
Mgiant = 1.00 + 0.07 Mo consistent with the earlier estimates ($4.3, 
Table 4). 

We also identify a significant blue veiling component, both through 
line veiling and second light contributions to the B and V-band light 
curves (§ 4.4, Figure 9). The veiling component contributes ~63% 
and ~24% of the total flux in the B-band and V-band respectively. 
The SED of the veiling component decays rapidly towards both bluer 
and redder wavelengths, strongly inconsistent with a stellar SED. 
Given that we do not see eclipses in the light curves, we infer that 
the veiling component has to be diffuse. 

We find no evidence for a luminous stellar companion and can rule 
out both single and binary main sequence companions based on the 


SED and limits on eclipses (Mgingle < 0.80 Mo, Mpinary < 1.51 Mo) 
from the light curves (§4.5). The SED and the absence of continuum 
eclipses imply that the companion mass must be dominated by a 
compact object even if it is a binary (§4.5, §5, Figure 10). 

Once the spectrum of the red giant is subtracted, we also find 
evidence of Balmer Ha and Hf emission ($4.6, Figure 11). The 
morphology of the Balmer emission lines is complicated and its 
origin is unclear. Even though we observe eclipses of the Balmer 
lines when the dark companion passes behind the giant, the velocity 
scales seem too low to be associated with an accretion disk. We 
also detect this system in the X-rays with Lx ~ 7.6 x 1029 ergs s7! 
(L/Leaa~107°) using Swift XRT data (§4.7). 

The simplest explanation for the dark companion is a single com- 
pact object, most likely a black hole, in the “mass gap”, making 
V723 Mon the host to the closest black hole yet discovered (§5, Ta- 
ble 5). Prior to this discovery, A0620-00 (V616 Mon) was the closest 
confirmed black hole at an estimated distance of ~1.6 kpc (Gelino 
et al. 2001). However, more exotic scenarios can also be plausible ex- 
planations, including a neutron star-neutron star binary and a white 
dwarf-neutron star binary. 

To better understand this unique system, further comprehensive 
multi-wavelength observations are necessary. In particular, UV ob- 
servations from the Hubble space telescope will constrain the nature 
of the veiling component, and X-ray light curves will be useful to 
understand the nature of the compact object in this system. Future 
data releases from Gaia will also confirm the orbital inclination. 

We can very crudely estimate the expected number of similar 
systems based on the fraction of the thin disk mass from which 
V723 Mon was selected. We assume a simple exponential disk model 
with density 


p ^ po e RíRa-czl/h. (14) 


where Ra « 3 kpc is the disk scale length (Amóres et al. 2017) 
and the numerical values of the disk scale height h and density 
normalization pg are not needed. The total disk mass is AnpohR?.. 
If we assume that V723 Mon was selected from a cylinder of radius 
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R at the Galactocentric radius of the Sun, Rọ ~ 8 kpc (Gravity 
Collaboration et al. 2019; Stanek & Garnavich 1998), the survey 
examined a mass of approximately 4z A R2 exp(—Rgun/Rq), so the 
fraction of the disk mass surveyed is approximately 


R 2 
r3 exp(- Rsun/ Ra) = 0.008 (15) 
Ra 


if we assume R ~ 1 kpc, as this encompasses most of the systems in 
the SB9 catalog. The fraction drops to 0.002 if we use the distance 
R = 0.5 kpc to V723 Mon. This implies that the Galaxy might contain 
some 100-1000 similar systems. Since ellipsoidal variability is only 
possible for a limited range of semimajor axes, it is not surprising 
that this estimate is larger than the number of mass transfer systems 
but smaller than estimates of the total number (10° to 10*) of non- 
interacting black hole binaries in the Galaxy based on population 
synthesis models (see, for e.g., Breivik et al. 2017; Shao & Li 2019). 

A number of large spectroscopy projects such as APOGEE (Ma- 
jewski et al. 2017) and LAMOST (Cui et al. 2012) are in the process of 
physically characterizing (kinematics, temperature, abundances, etc.) 
millions of Galactic stars. These surveys frequently obtain their spec- 
tra in multiple visits, providing sparse radial velocity (RV) curves for 
huge numbers of stars. A particularly important synergy is the ability 
to combine photometric surveys with these spectroscopic surveys to 
search for non-interacting compact object binaries like V723 Mon. 
For example, for the vast majority of these relatively bright stars, 
the ASAS-SN survey (Shappee et al. 2014; Kochanek et al. 2017; 
Jayasinghe et al. 2018, 2020) will supply all-sky, well-sampled light 
curves spanning multiple years. If we make a conservative assump- 
tion that ASAS-SN can characterize the variability of most giants 
up to ~3 kpc away (not accounting for extinction), there maybe 
~20 red giants with non-interacting companions that have ASAS-SN 
light curves. However, there is a significant cost to confirming these 
systems. In particular, a well sampled set of RV measurements is 
required to accurately measure the mass function and to constrain 
the properties of any companion. Nonetheless, as the spectroscopic 
surveys expand from a few 10° to a few 10 stars during the next 5 
years, this approach will become a major probe of compact object 
binaries. 
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Table 6. Swift UVM2 observations 


JD Date Phase UVM2 (mag) c (mag) 
2459144.439 2020-10-21 0.975 14.10 0.04 
2459150.607 2020-10-28 0.078 14.06 0.04 
2459155.040 2020-11-01 0.132 14.07 0.04 
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2459252.309 2021-02-06 0.774 14.03 0.04 
2459254.172 2021-02-08 0.806 14.04 0.04 
2459255.822 2021-02-10 0.833 14.06 0.04 
2459259.617 2021-02-14 0.896 14.05 0.04 
2459262.011 2021-02-16 0.936 14.03 0.04 
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Table 7. STELLA RV observations Table 7. STELLA RV observations, continued 

BJD Phase RV(kms!) ory (kms )) BJD Phase RV(kms!) cry (kms!) 
2454065.57926 0.242 —63.359 0.243 2455172.48258 0.709 64.109 0.091 
2454066.65035 0.260 —63.407 0.199 2455172.66102 0.712 64.566 0.118 
2454073.62965 0.376 —47.100 0.667 2455204.57655 0.244 —63.338 0.167 
2454073.66318 0.377 —43.862 0.305 2455208.54129 0.311 —58.753 0.153 
2454092.53351 0.692 61.584 0.199 2455209.59017 0.328 —55.705 0.195 
2454096.49062 0.758 66.656 0.104 2455210.52913 0.344 —52.216 0.203 
2454101.53952 0.842 58.195 0.166 2455211.48735 0.360 —48.608 0.205 
2454106.54614 0.926 34.825 0.288 2455212.58559 0.378 —43.778 0.195 
2454110.49315 0.991 8.326 0.633 2455213.43953 0.392 —39.561 0.214 
2454116.54293 0.092 —34.096 0.191 2455214.46909 0.409 —34.428 0.193 
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2454152.52001 0.693 61.471 0.142 2455231.52035 0.694 62.259 0.143 
2454162.44243 0.858 54.887 0.153 2455238.50895 0.811 63.532 0.172 
2454209.37594 0.641 50.398 0.692 2455239.46454 0.826 61.403 0.221 
2454359.72466 0.149 —50.362 0.294 2455240.50720 0.844 58.147 0.139 
2454364.72854 0.233 —63.249 0.217 2455266.45443 0.277 —62.580 0.148 
2454369.71403 0.316 —51.153 0.264 2455285.40569 0.593 36.069 0.148 
2454377.71952 0.450 —19.751 0.245 2455290.39565 0.676 58.852 0.114 
2454380.69831 0.499 0.271 0.21 2455297.3TTT2 0.793 65.281 0.144 
2454383.73052 0.550 20.050 0.215 
2454388.73313 0.633 48.436 0.228 
dae N bi iin Mie was made possible by the generous financial support of the W. M. 
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2454725.73812 0.256 —63.261 0.225 
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Table 8. Spectroscopic observations from HIRES, MODS and PEPSI. 


BJD [TDB] Date Phase 
2459143.09532 . 2020-10-20 0.952 
2459153.99899 . 2020-10-31 0.134 
2459162.02289 . 2020-11-08 0.268 
2459188.01112 2020-12-04 0.702 
2459208.01826 2020-12-24 0.036 
2459208.89679 2020-12-25 0.050 
2459209.98247 2020-12-26 0.068 


2459171.94982 
2459173.86921 


2020-11-18 0.434 
2020-11-20 0.466 


2459174.89034 2020-11-21 0.483 
2459175.77657 2020-11-22 0.498 
2459183.77896 2020-11-30 0.631 


RV (kms7!) pry (kms!) Instrument 
24.39 0.10 HIRES 
—47.28 0.10 HIRES 
—63.27 0.10 HIRES 
62.67 0.10 HIRES 
-12.23 0.10 HIRES 
-17.51 0.10 HIRES 
—24.92 0.10 HIRES 
— — MODS 

— — MODS 

— — MODS 

— — MODS 
48.74 0.25 PEPSI 
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APPENDIX A: THE ORIGIN OF THE RADIAL VELOCITY 
RESIDUALS 


$12 interpreted V723 Mon as a triple system where the companion 
is a binary with period Pi, = Porb/3 = 20 days. G14 argued against 
this hypothesis, but there are significant RV residuals at Porp/2 or 
Porb/3 depending on whether the model for the orbit of the giant is 
circular or elliptical (see $4.2 and Fig 3). We know from the limits 
on stellar companions in $4.5 that a star cannot dominate the mass of 
the companion, but there are also two independent arguments against 
any such binary, stellar or not. We do, however, have a hypothesis for 
the origin of these RV residuals. 

The first argument against a ~ 20 day inner binary is the dynamical 
stability argument raised by G14. For our nominal parameters, the 
outer orbit has a total mass of = 4Mg (in round numbers) and a 
semi-major axis of 100 Ro while the inner orbit has a mass of 3 Mo. 
The 45 Ro semi-major axis of a 20 day binary is very close to the 
Roche limit around the companion of 49 Ro. A rough estimate of 
the largest semi-major axis that could be dynamically stable given 
the outer orbit and the mass ratios is 32 Ro based on Mardling & 
Aarseth (2001). In short, G14 was correct to hypothesize that such 
an orbit should be dynamically unstable. 

We experimented with numerically integrating planar 3-body or- 
bits starting from nominally circular orbits as initial conditions for 
~ 100 orbits of the outer binary and varying mass ratios of the inner 
binary. High mass ratio inner binaries were generally very unstable, 
presumably because the lighter star in the inner binary is trying to 
maintain an orbital radius close to the full semi-major axis of 45 Ro 
in this limit. Equal mass inner binaries tended to be more stable, 
presumably because each star now only has an orbit of ~ 23 Ro 
about the center of mass of the inner binary. Nonetheless, many trials 
resulted in the destruction of the system well before 100 orbits were 
completed. These results were not unique to picking a truly resonant 
Porb/3 period for the inner binary. We suspect, but do not investigate 
here, that essentially all such inner binaries are unstable on long time 
scales. The inner binary would have to have a period shorter than 
— 12 days to satisfy the dynamical stability criterion, although such 
systems may still be secularly unstable. 

The other problem with longer period inner binaries is that they 
perturb the outer orbit and produce time varying tidal forces on the 
giant even if stable. In particular, the wider inner binaries would 
generally drive the outer orbit to be significantly elliptical in our 
numerical experiments even if they were stable over the 100 orbits. 
Making the companion a binary also means that the tidal forces on 
the giant are time variable. If the companion is a single star, the 
amplitude of the elliptical variability depends on (Rgiant/ a)? while 
if itis an equal mass binary it depends on 


Ro 2 1- R. 3 
x giant " x giant (Al) 
2\ di 2 d2 

where dı and d» are the distances to the two stars which comprise 
fractions x and 1 — x of the mass of the inner binary. Compared to 


a single companion, there is a fractional fluctuating tidal amplitude 
for circular orbits of order 


2 

; 15 a; 
3(1 = 2x) 9 cos wrt + = (1 - 3x +337) Lin cos2w,t — (A2) 

a 4 a2 
where the frequency wt = wj,—w is the frequency difference between 
the inner and outer binaries. In particular, for an equal mass (x — 1/2) 
inner binary, the peak-to-peak fractional change in the tidal forcing 
is 
15 "d 
2m. (A3) 
8 \a 
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so a Pin = 20 day period inner binary would produce ~ 38% peak- 
to-peak fluctuations in the tidal force with a period of 15 days. No 
such residuals are seen in the residuals from the ellipsoidal model 
of the light curves at this or any similar period (see Appendix B). 
To have fractional fluctuations in the tidal forcing smaller than f = 
0.1 fio requires an equal mass inner binary to be more compact than 
din < 23 fx, Ro and to have a period Pin < TAF days. 

These arguments appear to strongly rule out the inner binary pro- 
posed by S12. We instead suspect that the RV residuals and much of 
the evidence for a binary companion are driven by the consequences 
of making RV observations of a high amplitude ELL variable. This 
effect has previously been discussed (Hill et al. 1989; Eaton 2008) 
in the context of contact/semi-detached eclipsing binaries and el- 
lipsoidal variables in close orbits (Wilson & Sofia 1976). The very 
similar issues for observations of stars with dynamical tides have 
also been discussed (e.g. Arras et al. 2012, Penoyre & Stone 2019) 
recently. 

A tidally locked star is simply rotating at the orbital frequency, so 
the rotation velocity at the surface scales with the cylindrical radius 
from the rotation axis. Thus, the rotational velocity is larger on the 
long axis of the star than on the short axis. When viewed along 
a principal axis, there is the usual cancellation of the contributions 
from the parts of the star rotating towards and away from the observer. 
However, when viewed at an intermediate direction, the contribution 
of one sign of the rotation comes from the slower moving short axis 
while the contribution from the other sign comes from the faster 
moving long axis. This makes a contribution to the observed radial 
velocities with a period P,,,/2 because it is a m = 2 perturbation 
just like the tidal distortions. If the orbit is genuinely circular, a fit to 
the RV curve will instead yield a small eccentricity to try to remove 
the contribution from this effect (see Eaton 2008). This should mean 
that the dominant velocity residual will now have a period Pop /3, 
and we think this is the likely origin of the Po,5/3 signal found by 
Strassmeier et al. (2012). 


APPENDIX B: SHORT TIMESCALE VARIABILITY 


We searched for additional variability on three broad time scales. 
First, we looked for extra variability on the time scales of weeks or 
longer in the residuals of the ASAS and KELT light curves after sub- 
tracting the best model for the ellipsoidal variability. Next, we looked 
for variability on the time scale of days in the residuals of the TESS 
light curve. Finally, we searched for very short time scale variability 
during several periods of very high cadence ROAD observations. 

As discussed in Appendix A, a sufficiently wide binary secondary 
will produce strong fluctuations in the strength of the tides on the giant 
at periods with frequencies that are either the difference in frequency 
between the inner and outer binary w; or twice that frequency, with 
the perturbations at w; requiring an un-equal mass binary companion. 
In particular, for the 20 day period inner binary proposed by S12, 
we would expect strong fluctuating tides on periods of 30 and 15 
days, although the 30 day period could be partly absorbed into the 
ELL fit. We looked for periodic signals in the residuals of the fits to 
the ASAS and KELT light curves using Period04 (Lenz & Breger 
2005), and found none that were significant (a signal-to-noise ratio 
> 5). Given the level of the residuals and the expected amplitude of 
the perturbations for a 20 day inner binary, this appears to strongly 
rule out such a companion. 

We see clear evidence of structure in the TESS residuals on short 
time scales (Figure B1) with an amplitude of 0.4%. At least in this 
short TESS light curve, we find no evidence for a significant periodic 


signal on a time scale of days. These time scales are interesting be- 
cause the roughly correspond to the expected frequency of maximum 
power for asteroseismic oscillations. This frequency is defined by 
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where Ymax,o = 3100 wHz, Tego = 5777 K and go = 2.7 x 
10^ cm? /s (Brown et al. 1991; Kjeldsen & Bedding 1995). Given the 
spectroscopic parameters in $4.1, we estimate that vmax « 6.8 uHz, 
which corresponds to a period of Posce « 1.7 days. With no clear peaks 
in the periodogram of the existing TESS data, we cannot presently 
use asteroseismology as an additional probe of the mass of the giant, 
but this may change with further analysis and additional TESS data. 
Finally, many compact-object systems (e.g., Wagner et al. 1992; 
Zurita et al. 2003; Shahbaz et al. 2013) show rapid variability on time 
scales of minutes to hours. We used the ROAD telescope (Hamb- 
sch 2012) to continuously observe the system for about 2 hours on 
multiple nights, obtaining in total more than 3,000 10-sec B-band 
exposures. We have reduced these data using the ASAS-SN version 
of the Alard (2000) image subtraction software. We do not find any 
evidence for short («1 hour) timescale B-band variability, with the 
rms scatter during individual nights below 0.01 mag. Given that the 
"second light" contribution to the total B-band flux of V723 Mon is 
about 60%, this translates to the short-timescale variability of this 
component being less that ~ 2%. However, on longer timescales of 
hours to days B-band variability at the ~ 2% level is still possible 
(on top of the observed ellipsoidal modulations discussed in $4.3). 


This paper has been typeset from a TEX/IATEX file prepared by the author. 
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Figure B1. The light curve residuals for the TESS sector 6 data after the PHOEBE model is subtracted. 
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